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Abstract: We define matrix models that converge to the generating functions of a wide variety 
of loop models with fugacity taken in sets with an accumulation point. The latter can also be 
seen as moments of a non-commutative law on a subfactor planar algebra. We apply this 
construction to compute the generating functions of the Potts model on a random planar map. 

1. Introduction 

Loop models naturally appear in a variety of statistical models where the loops represent the 
configuration of boundaries of some random regions. In this article, we restrict ourselves to loop 
models on random planar maps. Perhaps the most famous of these is the so-called 0{n) loop 
model which can be described as follows. 

Consider the two-sphere 5^, and fix r disjoint disks 1 < j < r inside of S^. Each disk is 
given an even number of boundary points, one of which is marked. By a tangle we mean (the 
isotopy class of) any possible collection of non-intersecting strings in S"^ \ {'2i\ U • • • U 'Sir) joining 
the boundary points. 
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By a vertex (also called a Temperley-Lieb diagram) we mean an (isotopy class of) arbitrary 
non-intersecting collection of strings drawn inside a disk with a given even number of boundary 
points (one of which is marked). 

By a configuration P built on vertices Di, . . . , Dr wc mean a tangle Q into which the vertices 
Di,. . . , Dr have been inserted (so that Dj is inserted into in a way that boundary points 
match, and marked boundary points match). In the 0{n) model, the vertices are all assumed 
to be copies of the same vertex D: Let D be given by the following picture: 

The outer boundary is represented by a thin line. The boundary contains 4 boundary points, 
which are joined inside of the disk by non- intersecting strings (represented by thick lines). 
One of the boundary points is marked by a * to distinguish it from the others. A possible 
configuration is drawn below: 




Fix a number n (weight loop, sometimes also denoted S and called fugacity) and associate to 
a configuration P the value of P, computed as follows. For each , remove its outer boundary. 
We are left with a collection of closed loops (formed by strings outside and inside ^i, . . . , ^r)- 
The value of P is then given by 

value(P) =n# ^°°p" "P. 
In studying this loop model, one is interested in understanding the partition function 

oo 

/"(*)=Eh E Value(P) (1) 

r=0 ■ PeP(r) 

where the sum is taken over the set P{r) of all configurations built on r copies of D (note here 
that the vertices and boundary points are labeled, so configurations corresponding to different 
matchings of the labeled boundaries of vertices and the tangle give rise to different terms in 
the summation). 

One can consider a slightly more complicated situation, where exactly one of the vertices T 
(sometimes called the external face) is different from D. For example, T can be taken to be 

m 

We can then consider the observable 

oo r 

^t''{T):=Y,- Value(P) (2) 

r=0 " PePir,T) 



Loop models, random matrices and planar algebras 



3 



where P{r, T) is the configuration based on r copies of D and one copy of T. 

The loop model we described has been widely studied, in part because of its connection with 
the critical Potts model. In the case that n takes an integer value, the observable =5^" can be 
related to a random matrix model as follows. Let < i < n, he n M x M independent 

GUE matrices whose entries have variance M~^. For a vertex T, enumerate its boundary points 
starting with the marked one clockwise. Then consider the pair partition 

^T- i ~T i iff points of T are linked by a string inside T. 

For a Temperley-Lieb element T with 2p boundary points, define the polynomial Pt in n 
non-commutative indeterminates {Xi^ . . . , Xn) by 

Pt{Xi, . . . ,Xn) — ^ Xi^Xi^- ■ ■ Xi^^ . (3) 

l<ii ,...,22p<n 
il=ik if l~Tk 

Then, Wick calculus (see also its large M limit given by Voiculescu's asymptotic freeness |Voi91[ 
IVDN92j ) shows that 



lim E 

M->oo 



^Tr(PT(Ai,.-- ,A„)) 



where Tr is the standard unnormalized trace, Tr(A) — '^Au. Moreover, one can generalize 
this relation to negative real numbers t which are sufficiently close to zero by instead consid- 
ering random matrices chosen according to a certain Gibbs measure (rather than the Gaussian 
measure). For each M — 1,2..., consider the random choice of n matrices Ai, . . . , An of size 
M X M according to the probability measure 



d^^-''\A^,. . . , A„) = ^ exp |-MTr {y,^ A] - t( 



4^)2 \dA^...dAn 



where dA\ . . . dAn denotes the Lebesgue measure on the space of n-tuples of Hermitian matrices. 
This is the so-called 0{n) matrix model. The specific choice of the potential ~t(^J2^=i ^1) 
corresponds precisely to our way of representing by a polynomial the diagram D 



appearing in the definition of the loop model observable ,3^. 

The random matrix model then computes the observable: according to |BIPZ78| (see also 
|GMS06j for the derivation of the asymptotics of multi-matrix models) if i < is small enough 
for any polynomial P in the set of polynomial of n non-commutative variables 



lim 



^Tr(P(Ai,... ,A„)) 



dtiT''{A,,...,An)^-T:{P) 



where, if T is a Temperley-Lieb element, and Pt is given by ([3]) 

t?nPT) = ^t"(r). 

This representation of the 0{n) model as the asymptotics of a matrix model is useful in many 
respects. It was used to actually compute the partition function fn{t), see [DK88llKos89bKS921 
IEK95l[BE09| . 
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Finally, we mention that in the context of Voiculescu's free probability theory, the observables 
tr^ and !7 have an interpretation as a non- commutative law. Indeed, tr^ is a tracial linear 
functional on the set of polynomials in n non-commutative self-adjoint variables equipped with 
its natural product. For the observation is that polynomials of the form Pt form a subalgebra 
of the *-algebra of non-commutative polynomials in n self-adjoint variables. Multiplication of 
polynomials Pt, Ps corresponds to combining Temperley-Lieb diagrams in a disjoint fashion, 
TA S, and taking adjoints amounts to flipping the diagram along a vertical axis. The observable 
function =f7 then gives us a tracial linear functional on this algebra. 

Remarkably, for small real values of t, this linear functional is a state: tr^ [PP*] > for any 
P, a property which therefore also holds for 5^". This means that for "self-adjoint" diagrams 
T, the values of the observable {c$^"(T'^)}fc>o are moments of a probability measure, which 
turns out to be a useful observation. 

The non-commutative law r = tr^ can be called a free Gihhs state, the free-probability 
analog of a Gibbs measure. The combinatorics of this law can be compactly encoded by the 
so-called Scliwinger-Dyson (or loop) equation, which is a free analog of integration by parts, 
and for which the linear functional t — tr^ is the unique (suitably bounded) solution: 

t{XjP) = t (g) T{djP) - 1 t{p&jPd), yp (4) 

where dj and S>j are certain differential operators occurring in free probability theory |Voi98l 
IVoin2allW06ll(;MSn6l[Gmn9] . 

To summarize our discussion of the 0{n) model, we observe that we have (well-known) 
equality of the following values, for any TL diagram T: 

1. The value of the observable c$^"(T) obtained by summing over all possible configurations 
certain values associated to each configuration; it is defined for small values of t and all 
(possibly non-integer) n; 

2. The value of the free Gibbs law t{Pt)] here r was defined only for integer n as the unique 
trace on the algebra non-commutative polynomials in n variables which solves the Schwinger- 
Dyson equation (0]); 

3. The limit tr^ (Pt) as M oo of expected values E^„,m (M^^Tr(PT)); it is defined only for 
integer n averaging with respect to certain Gibbs measures on n-tuples of Af x M self-adjoint 
matrices. 

The goal of this paper is to generalize and extend these three objects by making use of 
the notion of a Jones ' subfactor planar algebra. The planar algebra is an abstract object that 
simultaneously encodes the value of the fugacity n and the possible choices of the vertices T, 
D in the loop model. 

Thus given a planar algebra V and having chosen elements S*, 5*1, . . . , 5„ of V, we construct: 

— A observable 5^ for the associated loop model; 

— A family of Gibbs measures /if^ defined on certain spaces of matrices with size going to 
infinity with M; 

— An algebra A = {Pt : T S V} with the linear functional ^ : A — >■ C satisfying a certain 
Schwinger-Dyson equation and a map A 9 T i— > Pt into our space of matrices. 

Moreover, we prove that these three objects are related by equalities: 

.%iS) = lim jj ( Tr{Ps)dtif'' . 

M^oo M J 

In other words, we introduce generalizations of loop models whose vertices are from an 
arbitrary planar algebra V; find analogs of the Schwinger-Dyson equation and of free Gibbs 
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states on graded algebras associated to V] and construct random matrix models that give the 
appropriate free Gibbs state in the large- Af limit. 

In particular, by specializing to the Temperley-Lieb planar algebra V = TL, our work 
allows one to represent the observable 5^" as random matrix integrals for non-integer values of 
n £ {2 cos(7r/p) : p > 3} U [2, -|-oo). This is important because -S^^ is analytic in n on a certain 
domain (see Lemma [5] in the appendix) and thus knowing it on a set with an accumulation 
point determines the function. Our work thus put on a firm mathematical ground the (often 
assumed) analytic extension of =3^" as a function of n, and this in a remarkably general setting. 

The passage to general planar algebras also naturally leads us to consider shaded tangles 
and vertices. This in particular allows us to consider a kind of U{n) model (closely related 
to the 0{n) model) in which configurations are shaded. This model corresponds to the Potts 
model, rather than just the critical Potts model described by the 0(n) model. Using our random 
matrix model, we are able to compute a observable of the U {n) model, giving a refinement of 
the computation in the 0(n) case. 

In the next three sections, we give a few details about our construction. 

1.1. Loop models associated to arbitrary Jones planar algebras. We begin by noting that it is 
possible to define an analog of the loop model if one permits the vertices T and D to lie in 
an arbitrary Jones subfactor planar algebra. We will review the definition of a planar algebra 
later in the paper; to give an informal idea of the connection we again consider a planar tangle, 
but this time also requiring that there be a checkerboard shading of the regions between the 
strings. The main axiom of a subfactor planar algebra then specifies that given such a tangle 
and k arbitrary elements of the planar algebra, one obtains a number, which is the value of 
the "configuration obtained by gluing in the vertices into the tangle" . Thus if we choose as our 
vertices an element T and n copies of an element D of the planar algebra, then one can make 
sense of the observable 5^ (T, D) just as before, by forming the summation over all tangles of 
the values of the resulting configurations. 

A particular example of a planar algebra is the (unshaded) Temperley-Lieb planar algebra 
TL{n), whose elements are exactly the Temperley-Lieb diagrams (the fugacity n is part of 
the planar algebra structure). Considering this planar algebra leads to the observable 5^" we 
described earlier in the introduction. 

Another example is the shaded Temperley-Lieb model, in which we require that the vertices 
be Temperley-Lieb diagrams with a fixed checkerboard shading of the regions between strings. 

Requiring tangles to be shaded comes naturally from the theory of subfactors: the planar 
algebra of a subfactor iV C Af is necessarily shaded, the shaded and unshaded regions corre- 
sponding to N and M respectively ^Jon99) . But planar algebras do not need to be shaded. In 
fact the original Temperley-Lieb equivalence jTH71| was a mathematical equivalence between 
a shaded model (the Potts model) and an unshaded one (Lieb's ice-type model). This equiv- 
alence is clearly described in chapter 12 of Bax82 where we see that the spins of the Potts 
model live in the shaded regions of a shaded 4-valent planar graph, while for the ice-type model 
there are only two "spin" states which live on the edges of that graph. A more general form of 
TL equivalence is that all representations of the Temperley-Lieb planar algebra will yield the 
same answer for the partition function of a model whose Boltzmann weights are defined using 
elements of the TL algebra. (The interactions may involve more than two spins, even in the 
shaded case.) 

In this paper we give matrix models for shaded planar algebras, for which 't Hooft diagrams 
in the perturbative expansion are only allowed if they can be shaded. In the particular case of 
the planar algebra underlying the Potts model, a quartic potential that uses the shading in 
an essential way is shown to yield the generating function for the Tutte polynomials of planar 
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graphs. It is conceivable that one may obtain the same generating function from an unshaded 
planar algebra by some kind of TL equivalence. But there are shaded planar algebras which 
are not equivalent to any unshaded planar algebra, and matrix models whose potentials are 
elements of these planar algebras would be inaccessible to unshaded models. 



1.2. Random matrix models. In the case of our specific loop model and for integer values of 
fugacity, the partition function of the model could be expressed as a limit of random matrix 

integrals. For example, if _D = ^ one is led to the random matrix model associated to 

the Gibbs measure on matrices given by ([T]). Considering a shaded analog of our loop model 
(corresponding to the observable given by 

arS (rp\ _ H ^2 r# loops m P /r-, 

2^ (ri!)(r2!) ^ 

ri,r-2=0 \ >- ^ PeP{ri,r2,T) 

where the sum is over all possible sum over all possible configurations involving a vertex T, ri 



copies of ^ ^ and r2 copies of 



) leads to a kind of U (n) matrix model which corresponds to 



integration over M x M non-Hermitian matrices , . . . , An with respect to the Gibbs measure 



1 I " 

d^ir,,t2 (Ai , . . . , A„) = ^ exp - M Tr ( ^ A, A* 

tl,t2 \^ j — l 



- J2 {tiAA*A,A* + t2A,A*A,A*) Ud^i . . . dA„. (6) 

i,J = l J 

It is not hard to generalize the definition of Pt given in ([3]) by replacing some of the Ai by A* 
to prove that 

aITc« (]^Tr(PT(Ai, . . . , A„))) = :^tlt2 (T) ■ 

It is natural to wonder if there are any random matrix models that can be defined for non- 
integer values of n. To be consistent with subfactor notation we will from now on use the letter 
5 to denote the fugacity (or loop weight) instead of n. 

In this paper, we show that this is the case for arbitrary loop models coming from arbitrary 
planar algebras (satisfying a technical restriction). In particular, this allows us to construct 
random matrix models that correspond to loop models with S in the set {2cos(7r/p) : p > 
3}U [2,-l-cx)). 

A key step in this construction uses the fact that arbitrary subfactor planar algebras can 
be viewed as subspaces of the linear spaces generated by loops on certain bipartite graphs (the 
so-called principal graphs). Thus in our matrix model the random matrices are indexed by the 
edges of a bipartite graph G ~ {V, E) with vertices V and edges E, whose adjacency matrix 
r has Perron-Frobenius eigenvalue 5, which restricts ourselves to the above possible values. 
As G is bipartite, V — U V~ . We will see that the configurations can be indexed by the 
vertices and the (oriented) edges of the graph; shaded (resp. unshaded) regions will be labelled 
by vertices from l/_ (resp. V+) and the strings between a region indexed by u and another by 
V will be indexed by edges between u and v. At the level of the matrix model, if A,, denotes 
the matrix indexed by the edge e from vertex s to vertex t in the bipartite graph, then A^. 
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will be a [/i(s)Af] x [/i(i)M] matrix with M a parameter going to infinity and ^J,{v),v £ V the 
Perron-Frobenius eigenvector of F for the eigenvalue S (here [•] denotes the integer part of a 
number) . We will require that A^o = A* where e° denotes the edge e with opposite orientation. 
To simplify, let us consider the planar algebra generated by shaded Temperley-Lieb elements. 
In that case the graph G can be chosen to be the Dynkin graph A,„ (with m a function of S) . 
Let B be a shaded Temperley-Lieb element. The generalization of the polynomial Pb of Q is 
then given by 

(JB{ei ■ ■ ■ e2k)Xei ■ ■ ■ ^e2fe 

ei=e° if J~p 

where we sum over all loops w = ei ■ ■ ■ e2k starting (and finishing) at u in G so that if j P, 
Cp is the opposite edge e° of Cj. Here CTB(ei • • • e2k) is an appropriate constant (see equation ©) 
which only depends on the Perron-Frobenius vector /i and the loop ei, • • • , e2k- Here v £ V-^ 
iff the marked point of B is in an unshaded region. 

To give a more precise idea of our matrix model let us first consider the case without 
interaction (that is when the parameter t vanishes) and denote = £^o,o with 5o,o defined in 
dS])- In this case, the entries of A^ will be independent Gaussian variables with covariance given 

by ^y ix{s)ii{t)M'^ . We denote /i^^ the law of the matrices Ae,e G E which are independent 
modulo the symmetry constraint Ago = A* . If we have a shaded Temperley-Lieb element S 
with marked point in an unshaded region, it follows from Wick calculus (see also |GJS07) ) that 
for any v £ V+ 

Hence ^o{S) can be constructed as a limit of matrix models. This fact was used in |GJS07) to 
construct tracial states on planar algebras. In the case with interaction, t ^ 0, we need to add 
an interacting potential as we did in the case where 5 = n is an integer number. This potential 
is defined in the same spirit as the polynomials Pg : namely, if we want to add in the interaction 
a shaded Temperley-Lieb element S with marked point in (say) an unshaded region we may 
consider the potential 

WsiX,, e£E)^ M(t')-Ps(^e, e £ E). 
Let us now define, for Si, . . . ,Sr shaded Temperley-Lieb elements the Gibbs measure 
diJ.fl_t^{A,,e £E)^ exp [ M Tr { ^ t^Wg, (A, e e S)} j d^l^' {A,, e £ E) 

tl,...,tr \ fe / 

(for this to be well defined, we may need to add a large enough cut off, see section [3.3.11 we will 
also require that the graph G be finite (otherwise we can show in this case that the measure 
can still be defined due to the exponentially small correlations of matrices with "far enough" 
edges, see section [2.6.2p ). Note that even though WsiA^^e £ E) is a, sum of square matrices 
with possibly different dimensions, its trace is well defined by linearity. 

For instance, the matrix model associated with the U{6) model ^ we already discussed is 
given, see ([S]), by 

rf<,t.(^) = ^exp<jMX;(i2Wy_+tiWy+)M«)Tr( ^ <y{e)A,A:\ \d^l''{A) 
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This Gibbs measure is well defined provided ti,t2 < 0. 

In the U (n) model for n integer, Temperley-Lieb diagrams can be represented as loops on the 
graph consisting of two vertices and n edges between them: to a diagram T with 2r boundary 
points we associate the sum of loops 

<^T = ^ CT(ei • • •e^jje^ieia ■ • • ei2,_iej2^ (7) 
% =e°, if j~Tfc 

where Cfc denotes the fc-th (oriented) edge from the first vertex to the second one, and 
denotes the same edge with opposite orientation. Thus the random matrix model involves non- 
self-adjoint M X M matrices Ak corresponding to the edge Cfc (the adjoint Al. corresponds 

to e^). The diagrams * ^ ^ and *|^^^ then correspond to the terms J^ij A*AiA*Aj and 

which (after a cyclic permutation not changing the value of the trace) can be 
recognized as the two quartic terms in the potential of the U{n) model. (Note that in the 
case where the fugacity takes integer values, because of the special structure of the graph, the 
Perron-Frobenius eigenvector is identically equal to 1 and so cr(ei • • • e^^^) = 1.) 

One of the main theorems of this paper is the following generalization of the representation 
of the observables associated to loop models by the asymptotics of random matrix models: 

Theorem 1. Let 6 € {2 cos(^)}p>3 U [2, oo[ and let S, Si, ... , Sk be k fixed shaded elements of 
the Temperley-Lieb planar algebra with fugacity 5. Forti, 1 < i < k, small enough real numbers, 
consider the observable 

k 

ri,...,rki=l ' PeP(ri,...,rfe) 

where P{ri, . . . , rfe) is the set of configurations build on S, ri copies of Si, r2 copies of S2, etc. 
and 

Value{P) = 5*'='°^'='^ 

Then, for all v (z V-^ (resp. V-) if the marked boundary segment is in an unshaded (resp. 
shaded) region. 



.%{S)^ lim 



^ Tr(P^(A„ee£;)) 



The same theorem holds for more general planar algebras (rather than the Temperley-Lieb 
one). This theorem is proved in section [21 



1.3. Graded algebras and free Gibbs states. We saw on page 2] that the observable ^ has an 
interpretation as a linear functional on the set of polynomials of the form Pt, where T is a 
Temperley-Lieb diagram. It is not hard to see that such polynomials form a subalgebra of the 
algebra of all non-commutative polynomials. In fact, Pt ■ Ps = Ptas if we define T A 5 as 
the disjoint union of T and S drawn so that their marked boundary points are at the top of 
their disks, and then inscribed into a larger disk. It turns out that the operation A still makes 
sense whenever T, S are elements of an arbitrary Jones planar algebra. The resulting associative 
algebra is denoted by Gro{'P). In consequence, the observable =^ can be viewed as a trace on 
GroiV). 
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Remarkably, we show that this trace is also positive-definite (i.e., it is a tracial state): ,%{T /\ 
T*) > for all T. Furthermore, the combinatorics of this trace is also governed by a kind of 
Schwinger-Dyson equation, which can be represented diagrammatically (see Lemma [1]). The 
differential operators appearing in this equation are analogs of the free difference quotient and 
the cyclic derivative from Voiculescu's free probability theory |Voi02b[IVoi02a[rVoi06| , and turn 
out to be important for subfactor theory (see e.g. [CJSllj ). 



1.^. Computations for certain models. In the last part of the paper, we turn to two classical loop 
models, and show that we can indeed use the matrix models we have constructed to compute 
partition functions of the loop models. Since in our random matrix models the fugacity 5 can 
take its value in the set {2cos(7r/p) : p > 3} U [2,+oo), this allows us to determine these 
partition functions for any fugacity by analyticity (see Lemma [5]). 

The first model we consider is constructed with Temperley-Lieb elements with non nested 
strings and black inside (that is, depending only on a cup shaded black inside with the notations 
of (GJSOTj) 




In section 13.11 see Lemma [31 we identify the law of cup (the element made with only one 
string and black inside) under as a probability measure minimizing a certain entropy, whose 
Cauchy-Stieltjes functional can be computed. 

The second model we consider is the one mentioned earlier on, built on two Temperley-Lieb 

elements with two strings and opposite shading 5*1 = ^ ^ and 5*2 =j^^^|. We study the law 
of cup under 3^: 



H H 

i' /' / ' 




By our construction, we can see this loop model as a limit of matrix models for any 6 G 
{2cos(7r/p),p > 3} U [2, +00). We show in section [3731 that this model is related to an auxiliary 
model which shows up thanks to the Hubbard-Stratonovitch transformation, see Proposition 
m This auxiliary model also allows to compute the generating function of cup under ,%. We 



10 



A. Guionnet, V. F. R. Jones, D. Shlyakhtcnko, P. Zinn-Justin 



then solve this auxihary model based on the remark that it depends only on the eigenvalues 
of the matrices involved. Therefore, standard large deviations techniques IBAG97| can be used 
and the asymptotics of this model are described by a variational formula, see Proposition [5] 
We finally study the optimizer of this variational formula and show that its Stieltjes transform 
can, up to a reparametrization, be expressed as a ratio of theta functions, see Proposition [S] 
We summarize below our main results on the Potts model. 

Theorem 2. Let S/'^ be the tracial state built on the two TL elements Si and S2 with two 
strings and opposite shading, Si having one unshaded region, see ([5]). Assume that t = (<i,i2) 
are non negative and small enough. 

— There exists two unique auxiliary probability measures v^) on the real line which min- 
imize 



+ \/8tix + \/8t2y\di^+{x)di'^ {y) 



with S the free entropy S{fi) ~ //log l^^ ~ y\dlJ-{x)d^{y). 

— Let Bn be the shaded Temperley-Lieb element with n non nested strings and black shading 
inside and put 

n>0 

The above series converges absolutely for 7 small enough. 

Let M^{z) — / X]n>o ^"^"'^'^+(^)' Then, ^*{z) — jzr^^rj^ is invertible on a neighborhood 
of the origin, with inverse z^^j), and we have on a neighborhood of the origin 

Ci,,t) = ^MHz\,)). 
7 

— The probability measures can also be described as follows. There exists numbers 
ai < 02 < bi < 62 so that i/^ (resp. v^_) is supported on [—02/ ^/8ti, — ai/V^ii] (resp. 
[{bi - (62 - l)/^/8^^]/ Moreover set 

dz' 



u{z) := - ai)(&2 - 02) 



62 \/{z' - ai){z' - a2){z' -~ bi){z' - 62) 



let z(u) be its inverse, and set 



Then, if q is such that 5 = q + q^^ , q ~ e*'^", we have 

uj+{u) = —^—j[{(p+{u) - (p+{-u)) + R{z{u))] 

with, if O is the theta function given by (|43p. 

, , 0(u - Woo - 2i^K) 0(u + Moo - 2i^K) 
(p+{u) = c+ — ^ +c_ — — ^ 
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and 

R{z) 




1 



The constants ai, 02, 61, 62, "oo, C-i-, c_, if are defined by implicit equations. If 6 — 2cos(^) 
for some p G N. ly is an integer and satisfies an algebraic equation. 

Whereas the matrix model for the first model is well-known, the matrix model for the second 
(a kind of U{n) model) has only been solved in some special cases in the literature, such as the 
0{n) model pK88,Kos89b,KS92,EK95,BE09 which corresponds to the case where the shading 
is neglected, that is ti = t2. Moreover, matrix models are usually provided for S integer, whereas 
our approach allows a general construction of the matrix model for all the values of S above, a 
set which accumulates at 2. For this specific model, one could a priori define a matrix model 
for all S by using the auxiliary model which appears thanks to the Hubbard-Stratanovitch 
transformation of section 13.31 this is however a rather indirect and specific approach whereas 
our matrix models apply in a much greater generality, in particular in situations where the 
model has not yet been computed. 

Our construction is closely related to Pasquier's jPasST) . though the latter is in a slightly 
different context, namely that of statistical models on the square lattice (whereas there is no 
underlying lattice in our construction; it is in some sense random). See also |Kos89ailKos92l 
IKos96| for another application of Pasquier's construction in the context of matrix models. 
Moreover, the enumeration problem corresponding to our second matrix model was recently 
considered in pBBMOQ in an equivalent formulation, namely the Potts model; we shall comment 
on the exact relation to our work in section 3. 



1.5. List of notations. 

S Loop weight (fugacity), or the loop parameter (part of the data of a planar algebra). To be 
consistent with established notation, we use S and n (see below) interchangeably f i32.1l p. 

[13. 

/ Partition function of a loop model. 

r = {V, E) A bipartite graph whose set of vertices is V and set of edges is E. The set V is the 
disjoint union of the sets V+ and V- of even and odd vertices f 8 32. 21 p. [14]). 

Gro{'P) The graded algebra associated to a planar algebra V. In the Temperley-Lieb case, the 
elements of Grg are Temperley-Lieb diagrams (drawn with all of their marked boundary 
points on the top, starting with the marked boundary point). The multiplication is then 
obtained by drawing diagrams next to each other and embedding the resulting diagram into 
a disk. (gOl p. E) 

Loops (closed paths) on a bi-partite graph; + and — indicates whether the loop starts at 

an even or an odd vertex ( H'Z.'Il p. ITU). 
IJ.{v) The value of the Perron-Frobenius eigenvector at a vertex v. ^2.21 p. [T4l) 
IJ,^ , /U^'^ /i^ is a Gibbs measure from which our random matrix model (of size proportional 

to M) is sampled. (glXH p. ^ 
M Matrices in our random matrix model have blocks with sizes proportional to M; we are 

concerned with the large M limit of the random matrix model. 
n The loop weight (fugacity) of a loop model when it is an integer number; the number of 

matrices or non-commutative variables. See 6. 
V Planar algebra. A reader unfamiliar with planar algebra may wish to concentrate on the 

Temperley-Lieb planar algebra TL, which can be thought of as the linear span of shaded 

Temperley-Lieb diagrams. ( ^12. 11 p. [T^ 
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Pk,Pk.± The fc-box space (fc-th graded component) of a planar algebra. The choice of sign 
± indicates shading. For the case of the Temperley-Lieb alagebra, Pk denotes all (shaded) 
Temperley-Lieb diagrams with 2k boundary points, and Pk,+ denotes those diagrams for 
which the marked boundary point is clockwise from an unshaded region. ft ^2.1[ p. I13p 
The planar algebra associated to a bipartite graph F. ( H'2.2l p. 

cr(e) The factor cr(e) — iJ.{t{e)) / ^{s{e)) associated to an edge e of a bipartite graph. 

aB For a planar algebra element B G V and an embedding i oiV into a graph planar algebra 
'P^ , o'siw) denotes the value of the image of i{B) on the loop w e . In the case of a 
Temperley-Lieb diagram B e TL, asiw) is given by dH) on p. [iTl 

S, T Elements of a planar algebra; particular examples of these are shaded Temperley-Lieb 
diagrams. 

s{e),t{e) For a oriented edge e, s(e) is the starting point and t{e) is the end-point. ( H2.21 p. 

1141) . We sometimes write s^w) for a loop w d L to indicate its starting vertex. 
t One or more parameters (e.g., ii, . . . , tk). 

Tt{A) The unnormalized trace of an L x L matrix A given by X^i^i -^a- 
Tto{S){v) For an arbitrary planar algebra 7^, the map ^{S) given by the pairing {S, X^TeTL 
between 5* and the sum of Temperley-Lieb diagrams gives rise to an element in . In the 
case of a graph planar algebra, ■S/'{S) is a function on the set of vertices, and Tro(5)(w) 
denotes the value of this function at a vertex v. i W2A[ p. [T5)) 
tr(A) The normalized trace of an i x i matrix A given by -^Tr(A). 

tr Linear functional obtained as limit of averages of traces with respect to a family Gibbs 
measures on matrices. 

T A non-commutative law, i.e., a tracial linear state on the algebra of non-commutative poly- 
nomials. 

An observable associated to a loop model; also a tracial linear functional on the algebra 

Gro{V) associated to a planar algebra V ( W2.'S[ p. [TE\i . 
TL The Temperley-Lieb planar algebra; the same symbol also sometimes refers to the set of 

all Temperley-Lieb diagrams, (fll p. ffiXH p. ^ 
V± Vertices of a bipartite graph, -I- indicates even and — indicates odd vertices. ( ^2.21 p. [H]) 
W The interaction term in our models; we assume that W = J^^i^i with Wi elements of a 

planar algebra V. 
II ■ I loo Operator norm. 

Acknowledgements. The authors are grateful to the referees of the paper for their helpful sug- 
gestions for improving the readability of the paper. 

2. From planar algebras to matrix models 

In this section, we introduce our matrix models and prove Theorem [T] The matrix models de- 
pend on a bipartite graph which shows up in Jones's subfactor theory (see [JPIO] and references 
therein). We first recall the definition of a subfactor planar algebra and the construction of the 
planar algebra from a bipartite graph. The example that the reader can keep in mind is the 
Temperley-Lieb algebra. 

We then define a family of random matrices associated to a bipartite graph whose adjacency 
matrix has eigenvalue S corresponding to some Perron-Frobenius vector /i. 

We next consider the case of a finite graph and introduce the Gibbs measure associated to 
and prove Theorem [T] for planar subalgebras of graph planar algebras. Finally, we extend 
our construction to certain infinite graphs. 
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2.1. Planar algebras. We first recall some generalities on planar algebras. The reader is referred 
to |Jon99) or |GJS07) for a more extensive introduction. 

Recall |Jon99) that a planar algebra is a collection of vector spaces V = {Pk} endowed with 
an action of planar tangles. 

A planar tangle is a drawing consisting of an output disk Dq and some number of input disks 
Di, . . . , Dk in the interior of (fc > 0). Each disk has an even number of marked boundary 
points. On each disk, one of the boundary segments is marked and called the initial segment. The 
boundary points are joined by strings drawn in the interior of Dq and outside all Di, . . . , Dk] 
in addition there may be some number of closed strings not connected to any of the Di's. All of 
the strings are non-crossing. Lastly, some of the regions between the strings are supposed to be 
shaded, so that each string lies between a shaded and an unshaded region. Planar tangles can 
be composed by gluing the output disk of one tangle into an input disk of another tangle so as 
to match up the initial segments. In doing so, one must ensure that the numbers of boundary 
points and the shadings match. 

The main axiom of a planar algebra is the existence, for each tangle T with disks Dq, . . . , Dk 
as above, so that Dj has 2bj boundary points, of a multilinear map Mt ■ P^^ x • • • x P^^ — > P^^' , 
where aj = + if the initial segment of Dj is adjacent to a white region, and aj = — otherwise. 
The maps Mt are supposed to be compatible with the operation of composition of tangles and 
invariant under isotopy. Moreover, the vector spaces {P^} are equipped with an involution 
compatible with Mt in the sense that MT{f*) — M^^T){f ° /)* for any orientation reversing 
diffeomorphism Lp. 

A subfactor planar algebra is a planar algebra so that dim(P|5^) = 1. As a consequence we 
can define a sesquilinear form on each P^ by 



{A,B) = 




where the outside region is shaded according to ±. We also require that ( , ) is positive definite 
and Mti = Mt2 where Ti and T2 are the following two 0-tangles: 




Once Po ± have been identified with the scalars there is a canonical scalar 6 associated with 
a subfactor planar algebra with the property that the multilinear map associated to any tangle 
containing a closed string is equal to 5 times the multilinear map of the same tangle with the 
closed string removed. By positivity of the scalar product, S has to be positive and in fact it 
is well-known that the possible values of S form the set {2 cos(7r/p) : p = 3, 4, 5, . . .} U [2, 00) 
[Jon83) . 

2.1.1. Example 1: Temperley-Lieb algebra TL. It is not hard to see that the Temperley-Lieb 
algebra is a planar algebra. Indeed, given a planar tangle T and some elements . . . , Bk S TL 
one can glue these elements into the input disks of T. Next, one can remove all closed strings 
by replacing each closed string by the factor 5. What results is another TL tangle, which is the 
result of applying the map Mt to Bi, . . . , Bk- Clearly, Mt is invariant by isotopy and P^ has 
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dimension one. Finally, the canonical scalar product is positive definite according to |Jon99) . 
One way to prove it is by verifying that the map of TL into a graph planar algebra is a planar 
algebra map. It thus takes the canonical bilinear form on TL to the canonical bilinear form 
on a graph planar algebra, where non-negativity can be verified directly. We will write TL((5) 
when we wish to emphasize the loop parameter (fugacity) S. 

2.1.2. Example 2: The planar algebra of two stitched Temperley-Lieb algebras. Later in the 
paper, we will use (in addition to the Temperley-Lieb algebra) the stitched planar algebra 
V ~ TL((5i) © TL((52). This will be needed in order to realize the so-called 0{n,m) model in 
physics. 

The n-th graded component of P„ of V is given by 

where the sum is over all partitions tt of {1,2,..., 2n} into two subsets of even sizes 2(p7r) and 
2(q'^) and (TL © TL) J = TL^^ (g) TL^^ . Graphically, one can view P„ as the span of the collec- 
tion of isotopy classes of tangles obtained by superimposing an arbitrary TL tangle colored red 
with an arbitrary TL tangle colored black, in such a way that the strings of the red and the 
black TL tangles are allowed to only intersect transversally, and so that the resulting tangle 
has a total of 2n boundary points (counted regardless of color). The partition tt then corre- 
sponds to the colorings of the boundary points of the resulting diagram. We assume that the 
checkerboard coloring of the two superimposed TL diagrams are retained and are independently 
superimposed. 

The isotopies need not preserve the intersections of the red and black strings, but must 
preserve the partition tt. We also assume that one of the boundary regions is marked "first" 
(it could be of either color). Two different isotopy classes of diagrams, T and S, are presented 
below (black strings are indicated by solid lines and red by dotted lines, and the four possible 
shadings of regions are indicated by I I , ^H, EiSiSl, ^H) : 




The planar algebra structure of P„ is defined as follows. Given a planar tangle T and diagrams 
Ai, . . . ,Ak 'mV, the result Mt{Ai, . . . , Ak) is obtained by gluing the diagrams Ai, . . . ,Ak into 
the input disks of T and summing over all possible ways of extending the colorings and shadings 
of the Ai^s to the resulting tangle. The construction of is a particular case of a more general 
construction Vi © 7^2 , which is possible for any pair of planar algebras Pi, 1^2- This construction 
is presented in Appendix [U] 

2.2. On the planar algebra of a graph. Following |GJS07| . we shall use the construction of 
planar algebras from bipartite graphs, as introduced in | JonOO] . The key ingredient here is the 
fact that every subfactor planar algebra (in particular, TL) embeds into a graph planar algebra. 
This makes it possible to "coordinatize" planar algebras. 
We first fix notations. 

Let r — (y, E) be a bipartite graph with vertices V — V- U V-^- so that any edge is either 
from V+ to V- or V- to V+. We denote by E+ (resp. by E-) the set of (oriented) edges starting 
in V+ (resp. V-). Thus E = E+ U E^. We let /i be a fixed Perron-Frobenius eigenvector with 
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eigenvalue 6 for the adjacency matrix of F. The vector //, has positive entries. If e G -E, wc 
denote by e° the edge with opposite orientation. We denote by L the set of loops on 7^, L+ 
(resp. L~) the set of loops starting in V+ (rcsp. V-) (so L = i+ \JL~). We denote by L{v) the 
set of loops starting at v E V. We finally let = U„^±P,f _|_ where P^± is the vector space of 
bounded functions on loops on F of length 2n starting and ending in V+ for the plus sign and 
V- for the minus sign. In the following, s(e) (resp. t{e)) is the starting (resp. target) vertex of 
an edge e (note that several edges can have the same starting and ending points). 

Example 1. Consider the graph with one vertex in V+, one vertex in V- and n edges between 
them. In this case, (5 = n is an integer, and the eigenvector ji is identically equal to 1. The TL 
algebra embeds into the planar algebra of this graph for integer fugacity 5. 

Wc next describe the action of planar tangles on V^. Let F be a pi anar tangle with k input 
disks and let Li, . . . , be loops on F. To define the planar algebra structure, we must exhibit 
Mt{L\, . . . , Lfe) as an element of the planar algebra, i.e., we must prescribe its value on a loop 
L. The value 

MT{L^,...,Lk){L) 

is computed as follows. First, label the marked points on the input disks of T with the edges com- 
prising the loops Li, . . . , Lfe, clockwise starting from the marked vertex (and the beginning of 
Lj). Next, label the marked points on the output disk of T with the edges of L, clockwise starting 
from the marked vertex. As a result, we obtain a labeled tangle, and we'll set Mt{L\, . . . , Fk){F) 
equal to the value of this labeled tangle which we compute as follows. First, we remove all closed 
loops in the tangle T (let us say, p loops total) and multiply Mt{Li, . . . , Lk){F) by 5^. We'll 
again denote the tangle with removed interior loops by F. 

Next, we isotope the tangle in such a way that each input disk becomes a rectangle whose 
top is horizontal, and so that all strings emanate from the top (in this way, the marked initial 
segment comprises the sides and bottom of the rectangle). Put 

V M(s(e)) 

Then the value Mt{Li, . . . , Lk){L) is zero unless each string connects points which are labeled 
by opposite edges. Otherwise, 

MT(Li,...,Lfe)(L)= [] a(e,)"^=/- (8) 

strings s 

where eg is the start of s and Og = dO is the total winding angle of the string s. Here d9 
stands for the 1-form ydx — xdy on the coordinate plane. (Note that the choice of which edge 

is selected as the start of a string s is irrelevant: if s' is the string s traversed backwards, then 
we get a non-zero value for the tangle iff eg/ = e°; note that cr(e°) = cr(e)~^ and Ogi = —9s-) 

We note the identities (thick lines indicate an arbitrary number of parallel strings, e and / 
are two arbitrary edges, and x,y,z,w arc arbitrary paths on graphs so that xey and zfw are 
loops, and all planar algebra elements are arranged so that the marked initial segment is at 
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top- left): 



LI 



l_l 



f 


= <5e=/°(T(e) 












U_l 



X w z 



3 



where a{z) = Yl'^i^j) if 2 = zi . . . z;. (The last operation can be though of as replacing the 
string connecting e and / by a "tunnel" joining the two planar algebra elements: 



z w 



followed by an isotopy "straightening" the resulting picture). 

Using these operations, any labeled tangle can be simplified leaving only strings connecting 
the outer disk to inner disks. Each of these remaining strings can be removed by contributing 
a certain multiplicative factor according to ([5]). 

As mentioned in the beginning of the section, we extensively make use of the following fact 
(see e.g. |JonOO] ): 

Proposition 1. Let V be a subfactor planar algebra. Then there exists a bipartite graph F and 
a planar algebra embedding i : V . 

Note that if V is embedded into T'^, then any B ^ V can be regarded as a function on the 
space of loops of F. We denote by ct_b(w) the value of this function on a loop w & L. Naturally, 
asiw) depends on the specific embedding of V (indeed, a fixed planar algebra V may admit 
embeddings into graph planar algebras for many different F). 



2.2. L Example 1 continued: Temperley-Lieb algebra inside the planar algebra of a bipartite 
graph. A particular case of the planar algebra axiom is the natural embedding from TL(A:,±) 
into a linear span of loops, following |GJS07] . Indeed, Temperley-Lieb tangles are tangles with 
no input disks, and thus produce elements in any planar algebra. 

Suppose that we are given a box B with 2k boundary points (arranged so that all boundary 
points are at the top and * is at position from the top- left). Assume also that there are 
k non-crossing curves inside B which connect pairs of boundary points together. Let tt be 
the associated non-crossing pairing. We let Lb be the set of loops in 7^ so that w E Lb iS 
w = ei ■ ■ ■ e2k with 
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Sn = £° if {n, £} is a block of the partition tt (which is denoted n 
s(ei) e V+ (resp. in V^) if B G TL(fc, +) (resp. B e TL(fc, -)). 




For e e E, a{e) := 




and for e L, we define the weight 



<Jb{w) = cr(eij • • • cr(ei,J if Ci^ = whenever w £ Lb, ik ^ jk and ife < jk, (9) 
(Jb{w) = ii w ^ Lb- 

We then associate to the Temperley-Lieb tangle B the element G given by 



Informally, wb — X^^eL ''^b(^)^' where we identify a loop w L with the delta function 
supported on that loop (viewed as an element in V^). Of course, this is only correct when the 
graph is finite, so that the expression for is a finite sum. 

2.2.2. Example 2 continued: The planar algebra of two stitched Temperley-Lieb algebras realized 
inside a graph planar algebra. Let V = TL((5i) © TL((52) as in H2.1.2I We'll realize V inside a 
graph planar algebra. Let ErjEh be two graphs so that the associated planar algebras contain 
TL{dr) and TL{Sb), respectively. We can thus choose to be Aoo if Sx > "2 and otherwise A„ 
for some n (related to S^). By appendix [Ul V embeds into the graph planar algebra of F. 

Let r = Fr X Ff). More precisely, the vertices of F are pairs (u^ , Vb) with Vx S Fx, x G {b, r}. 
The pair {vr,Vb) is even iff either both Vr,Vb are even or both are odd; the pair {vr,vi,) is odd 
otherwise. The edges of F are of two kinds: the red edges, consisting of pairs (e, v) with e an edge 
in Fr and v a vertex in this is an edge from {s{e),v) to {t{e),v); and black edges, consisting 
of pairs {f,w) with / an edge in Fi, and w a vertex in Fr; this edge goes from {w,s{f)) to 
{w, t{f)). For an edge in I^, e = (/, w), put u{e) — f (which is in Fr or Ft according to whether 
e is red or black) . Note that is a bi-partite graph, since each edge in F changes the parity of 
one of the components of a vertex {vr, Vb). By Appendix [Cl F is the principal graph of V. 

Let /i be the Perron-Frobenius eigenvector for 7^, given at a vertex (v,w) by the product of 
the eigenvectors of Fr and Fb. For e and edge of F, put a{e) = (/i(t(e))/^(s(e)))^/^. Let c{x) 
be the color of the a;-th boundary point of T. 

Let r € Pri be a diagram, and let 



The embedding of V into the graph planar algebra of F is given by sending T to the function 
/t G given by: 



wb{L) = (Tb{L). 



(10) 



Rt — {ix,y) : boundary points x and y are connected in T}. 



friei ■ ■ ■ e2n) 



n ^(^)'^' 



e has same co] 



)lor as c(x) ^u(ej:)—u(ey)° , 



{x,v)&Rt 
x<y 



for any loop ei • • • 62,1 in F . 
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2.3. The observables S/'t. From now on, we shall fix a subfactor planar algebra V as well as an 
embedding of V into the planar algebra associated to a bipartite graph 7^. Moreover, we 
shall fix elements Bi, . . . , Bk G V . 

For ti,l < i < k, small enough real numbers and any S G V , consider the observable (see 
the discussion in ijl.ip 

E E A^T(5,i?i,...,Si,...,i?fe,...,B0 

ri,...,rki=l TeT(ri,...,rfc) ' ^I'J ' ' ' 

where T(ri, . . . , r^) is the set of all possible planar tangles (configurations) with 1 + ri + • • • + rfc 
input disks and Mt{S, . . . , Bi, . . . , Bk, . ■ . , Bk) denotes the value (in Vo) of the tangle T 
applied to S and ri copies of r2 copies of i?2, etc. 

Since we have assumed that is a subfactor planar algebra, T^g — C and so the value of 
5f(S') is a complex number. However, the expression for ^t{S) makes sense for any S S . In 
this case, Mt{S, Bi, . . . , Bi, . . . , Bk, ■ . ■ , Bk) is valued in Vq , i.e., it is a function on the set of 
vertices of F (this function is constant if S actually happens to be in T'^). 



2.4- Random matrices associated with a graph. In the sequel, we fix a graph F with an eigen- 
value S and a Perron-Frobenius eigenvector /i as in the previous part. For e G E^, is a 
[Ms(g')] X [Af((g)] matrix with i.i.d entries with variance (M^^g^Afj^g-))^ 2 for some integer num- 
bers {My ,v GV) so that 

iim — — = fi(v). 

We put A^o — [Ae)* and denote, for a word w — ei - ■ ■ Ck, A^j — A^^ ■ ■ ■ A^^. /i*^ is the law 
oi A = {Ae,e e E) chosen independently except for the last constraint. In order that A^^ is a 
square matrix, we shall assume that w is a loop. 

In the particular case that w is the trivial path (consisting of no edges) starting and ending 
at V, we will denote by Ay^ the identity matrix of size My x My. 

The center-valued trace Tro on Gr^V^ is given by the equation, for x ~ ei - ■ ■ e2k 

Tro{x){v) = ls{ei)=v {x,Tk) 

with Tk — X]_BgTL(A;) sum of all TL diagrams with k strings and where for two loops 

X, y on F, {x, y) — dx=y. Thus Tro(a;) is a complex-valued function on the set of vertices of the 
graph. 

We have the following theorem from |GJS07I Proposition 2] 
Theorem 3. Let v G V and w = 6162 • • • € L{v). Then 

Tto{w){v)= lim [ iv{Ay,)dti^' {A) (11) 

A'f— >oo J 

where tr is the normalized trace on matrices, tr(i?) = ^^^Yl!i=i^ii if R is a N x N matrix. 
In the case where we extend linearly Tro o.'^d take w = wt for some subfactor planar algebra 
element T, Tro^wr) — C{T)1 is constant and C(T) = 3q{T). 

We used above the notations of ([7]) and (|10p . Note that in |GJS07| . we had an additional 
dimension N so that X^'^ converges to free variables entries as N goes to infinity. This is 
however not needed since M goes to infinity. 

We shall give a new proof of Theorem |3] based on the so-called Schwinger-Dyson equation 
and obtain its generalization to the case of non-trivial interaction (t ^ 0): 
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Theorem 4. Assume that the subfactor planar algebra is the Temperley-Lieb algebra, or is a 
planar subalgebra of a graph planar algebra of a finite graph ( this is the case, for example, if it 
is finite- depth). For t = {ti, - ■ ■ ,tk) S K'^ small enough, there exists a probability measure jif^ 
on the set of matrices indexed by e € E with dimensions M^^g) x Mt(g-) so that for v Cz V and 
w = 6162 • • • Cfe e there exists a limit 

= lim / tr(A„)d/xf (A). (12) 

Af — foo J 

where denotes the monomial = X^-^ ■ ■ ■ X^y, ■ In the case where w = wt for some 
Temperley-Lieb diagram T (or, more generally, T comes from a subfactor planar algebra inside 
V^), with = 'Yl,weL(v)'^T{w)Xy^ , trt(P^)(w) = Ct{T) is constant for all w € V+ (resp. V^) 
if the marked point of T is in an unshaded (resp. shaded) region. Moreover, Ct{T) — S^t{T) is 
given by ([2]). 

In the statement above, trj is extended by linearity and we used the notations of ([7]) and (jlOp . 
In the next subsections we shall prove this theorem. 



2.5. The case of finite graphs. For 5 — 2cos(^), p > 3 (or, more generally, if the planar algebra 
under consideration is finite-depth), the graph T can be chosen to be finite. This is the case 
for TL((5) a 5 — 2cos(;j^); the graph is then A„, the linear chain of n vertices and n — 1 
edges. Another such example is the case of the graph with two vertices joined by n edges (see 
e.g. [JonOOl Examples 4.1 and 4.4]). We study the finite-graph case first. 

2.5.1. Definition of the matrix models. Let V C be a subfactor planar algebra realized inside 
a graph planar algebra. We consider the law /j,*^ of \E+\ independent Mg(e) x Mf^^^ matrices 

with i.i.d Gaussian entries with variance (Ms(e)Mt(e))~^. We denote by ||-R||oo the spectral 
norm of a matrix R (that is, the spectral radius of V RR*), and we consider A ~ (^e)ee e+ , the 
collection of these matrices and ||A||oo := maxgg^;^ ||Ae||oo- Let us assume that V C V'" and 
Bi, . . . , Bk G V are fixed. To fix notations, let us write each Bj as a linear combination of the 
basis for V^, consisting of all loops in the graph T: 

Bj = ffB, iw)w. 
weL 

(Note that liV = TL, then (w) is exactly given by the formula ©). 
For given real numbers ti, . . . ,tk and some fsT > 2, we set 

f,f'^{dA) -.^h^^expLlY^U J2 KsH)<yBAw)TTiAA /^(dA), 

where we denote by s{w) the starting vertex of a loop w d L. 

Example 2. Consider the TL algebra for 5 e {2cos(^),p > 3}. Then TL can be embedded into 

the graph planar algebra where E = A„. We can consider the following potentials: (i) Let 
B be the element of TL given by the tangle with only one string: 
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We can either see it as a cup with black inside and white outside or with the opposite shading, 
both leading to the same potential in the matrix model : 



Y Kv) E '^(e)Tr(A,A:) = Tr 

i^G V+ e:s{e)—v 



vev+ 

\s{e)—v 



J 



v^V~ e:s{e)—v 



(ii) Let B be the element of TL given by a tangle with two strings, two white regions and one 
black: 




Then, it contributes to the potential by 



E M(«)Tr 
vev+ 



E E <r{e)cy{f)A.AfA)Al 



s{e)—v 



f 



E Ai(«)Tr 



E <^)AeA: 

e 



(13) 

where it is understood that products which make no sense give no contribution. Inverting the 
shading amounts to exchanging V-^ and 



2.5.2. The main result. To state our main result, as before, we fix a planar subalgebra V of 
a graph planar algebra and, for some Bi G V, let W = J2 tiBi, ti G K. We let ^ be the 
observable functional considered in 



Proposition 2. (a) Let K > 2. Then there exists t{K) > so that whenever maxi<i<fc \ti\ < 
t{K), for any v G V and for any loop w G L{v) there exists a limit 



(b) The Schwinger-Dyson equation (described below in ^2.5.4\ l has a unique solution Tt. 

(c) For any element w € , Tt{w) — trt{Au]) = S^t{w). 

(d) Let V C V'" be a subfaetor planar subalgebra and assume that W G P. Then for any 
Q = Hw'^wW e V and any v e F, Tt{Q){v) = trt(X]„ = ,9t{Q){v) are constant 
independent of v. 

The convergence of the matrix model is a small generalization of |GMS06[ Theorem 3.5] (where 
only Hermitian random matrices where considered), whereas the identification of the limit 
is based on the analysis of the so-called Schwinger-Dyson (or loop) equations. Note that in 
most papers in the physics literature, the cutoff K < oo \s not taken, leading sometimes to 
diverging integrals. The advantage of adding this cutoff is that all integrals are well defined and 
moreover for small ii's, the Gibbs measure nf^'^ has a strictly log-concave density, providing 
many interesting properties which allow to put on a firm mathematical ground the above 
convergence. In fact, we can remove the cutoff in case the density is strictly log-concave. K 
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has to be chosen strictly greater than 2 (which is greater to the edge of the support of the 
Pastur-Marchenko law) so that the limit does not depend on it for t small enough. We start 
by recalling these properties and sketching the proof of ^GMS06., Theorem 3.5]. 

2.5.3. Convexity assumption and consequences. For a loop w € L, let us put s{'w) for the 
starting vertex of the loop. As in ^GMS06_, we shall restrict ourselves to the case where ' 
has a strictly log-concave density, which amounts to assume that the map from the set of \E^-\ 
Hermitian matrices into K given by 

k 

iAe)eeE+ ^ - E E (HTr(A^) + ^ E Ksie)Me)TTiA,Al) (14) 

is strictly convex (with Hessian bounded below by c for some c > independent of M) on 
{ll^lloo < K}. This is always true for t sufficiently small, depending on K, with c going to 
m — min{^(s(e))cr(e) — ^{t{e))fi{s{e)), e G E+} > as t goes to zero 

As a consequence of strict convexity, \ 
|Gui09[ Section 6.3] namely for any w £ L 



As a consequence of strict convexity, we have concentration inequalities under nf^'^ , see 



M,K 
Mi 



The family 



tr(A.)-Mf^(tr(A.))p)<^^^. (15) 
{/if^(tr(A„)),u;eL} 

is tight. We will denote by {trf"(X^) : w G L} a. limit point. We also have Brascamp-Lieb 
inequalities, see |Gui091 Section 6.5], and so by comparison to the Gaussian law for which we 
know that the spectral radius is bounded with overwhelming probability, we can prove that 
there exists £{c) (which only depends on c) so that 



S(c)M 



-K 



for some 5{c) > 0. We assume that we have chosen K > i{c) so that the limit point tr^ will 
not depend on K. We denote in short 

trt(eie2 • • • efc) = trf (Xe^Xg^ • • ■ J 

for any path 6162 • • • in G. By definition the above vanishes if s(ei) 7^ t{ek). We next show 
any limit point satisfies the so-called Schwinger-Dyson equation and that this equation has a 
unique solution when the i^'s are small. 

2.5.4. Schwinger-Dyson (or loop) equations. Let us fix e G ii' and P ^ u with u a path from 
t{e) to s(e). By using integration by parts and concentration inequalities, we obtain (see |Gui09l 
Section 8.1] or [GMS06t Theorem 3.1]), that 



lim 

M^oc 



M.K 



An ( M(s(e))cr(e)yle ^ E^* E t^(^M)"'Bi {w){Dx^oA^) 




i—1 wGL 



Tr\ M.K /Tr 



0, 
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where dx^o (resp. Dx^o) is the non-commutative (resp. cyche) derivative with respect to Xe". 
As before, we write s{w) for the starting point of the loop w. 

Recah that if a; is a path starting at a vertex w, then tr(a;) is a hmit point of 
seen as a hnear form on GroT'^. Using this we see that tr satisfies the foUowing Schwinger- 
Dyson equation: for every path u from t{e) to s(e), 

^(i(e))tr I u { ^{s{e))cr{e) e - ^ ^ ^{s{w))(Tb, {w)Deo (w) 

\ ^ i w 

~ Ai(s(e))Ai(t(e))F (g) E (deou) = (16) 

where 

w—wie°W2 w—wie°W2 

Moreover, note that as w G L, Dgow is a hnear combination of paths starting at s(e) and 
finishing at t{e). Finahy, considering the trivial loop ly at a vertex v, we have that tr(l^) = 1 
by construction. 

Dividing (jl6p by /i(s(e))/i(t(e)) and then by cr(e), we obtain the following equivalent equa- 
tion: 



) 



tr(we) — tr I 7/ 1 U a{e) 

\ ^ i w 



iM(£H), 
/i(s(e)) 



(TB,iw)Deo{w) ] +(T(e)-Hr(g)tr(5eou) (17) 



We can use tr to define a collection of linear maps from with values in , which we 
will for now denote by tvij. 

By definition, if P G V'" is the delta function on the loop ai, 02, . . . , a^, ei, . . . , ej, 61, . . . , 6^, 
then 

tri,j(f ) = <5s(6i)=t(a.)ti'(ei • ■ -ejOQ 

where Q is the delta function supported on the loop oi ... 0^61 ... &r- 

Then the Schwinger-Dyson equation is an equation on these linear maps. We will use the 
following graphical notation for the result of tr^ j- applied to the delta function supported at 
the loop ai, . . . , Oi, ei, • • • , e^, 61, . . . , &j, € (we suppress i,j since they can be read off from 
the numbers of the various strings): 



tr 



(18) 





ai 


ei 62 Cj 


fcl 








A 

















We have the following lemma. 
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Lemma 1. The Schwinger-Dyson equation (|16p is equivalent to the following diagrammatic 
equation, which is required to hold whenever some arbitrary element of is inserted into the 
blank boxes (same element to each blank box) and the strings at the bottom of each of the three 
pictures are connected to an arbitrary element of (same element for each picture): 



tr 




i odd 



where = Ei ^iBi Ei Eu; ^s. (w) 



E 



tr 



z 



© 
W 



(19) 



Here we use the following conventions. Thick lines indicate an arbitrary number of strings, 
e we'll set 



Also, given 



" ■ " 


_ 


X 





where by convention there is 



an arbitrary fixed number (same in all diagrams) of strings between the marked initial segment 
of X and the last string that is bent around x on the left. 



Proof. Consider e, /, ei, . . . , e„ G E{r) so that the paths ei, 62, . . . , e^, 17 and e, / form loops. 
Consider the following planar operation applied to P = ei, . . . , Ck, g and e, / (more precisely, 
P is the element of which is the delta function on the loop ei, . . . , ek,g, etc.): 

















e 


















(20) 


■9 




/ 
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Let ei . . . Ck g a„ . . . ai and e / be two loops in F . Consider the following equation: 




tr 



X 



tr 



z 





ei---ei_i e; 




E 






i odd 







tr 



E 







^ 




o 




























) 





(21) 



It is easily seen (by noticing that g and / are arbitrary subject to t{f) — s{g)) that (PT|) is 
equivalent to ((TO)) . 

Next, note that both sides of (PT|) are zero unless the region containing the point at infinity 
in all diagrams is labeled hy v = t{e) (so in particular t{ai) = s(a„) = v). Consider now the 
three diagrams comprising equation (1211) . Put 

R = Sf^goa{f)ai ...an- 

The diagram in the upper-left corner is exactly 

Qi = tr(ei . . . efee)i?. 

Consider now the diagram in the upper-right corner. The term in the summation is zero unless 
= e°. In particular, the region to the right of the string emanating from Ci can be labeled 
by s(e) and the region to the right of the string emanating from e^-i can be labeled t{e). By 
definition of tr.y , we then get that this diagram is equal to 



(32 = tr tr 



ei . . . ei-i ® Ci+i . 

\i odd 
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Note that we can actually replace the sum over odd i by the sum over all i, since tr(ei . . . ei_i) = 
unless i is odd. 

Assume first that W is the delta function supported on the loop gi . . . gs- Let us now consider 
for i even the diagram 



D, = 



W 




Since all strings emanating from gi+i, • ■ • ,5s make a 360° rotation in the drawing, according to 
(jS]) we have 

D,^<j{g,+,)^---a{gs)^Dr 

where 

^ 1 

D' = 



gi+i, 



To see this, we can compare the results of glueing all strings of Di to some diagrams . . . , 
versus the results of glueing the corresponding strings of D[ to the same diagrams. As we apply 
([Sj to remove strings, the strings starting from Di associated to gj,i-\-l<j<s contribute an 
extra factor of (j{gj)^ as compared to their contribution if they were to start from D'^. 
Using this, the value of the bottom diagram in (l?T|) is 



Qs = ^ cr(gj+i)^---cr(.gs)^tr(ei...efc5i+i...5sgi....gi_i)5e°=g,cr(e) 

i even 
i even 

m(s(.9i)) 



/i(s(e)) 



tTiei...ekDeoW)aiey'R 



By linearity, the same equation holds for arbitrary W: 

Qs = o-g. ("w) ^^'^l^NN tr(ei . . .ekDe'>w)a{ey^R 

^ ^ li\s{e)) 

% W r \ \ / / 

Now, ([21]) is equivalent to saying that 

Qi = Q2 + Q3- 

Since e, /, oi, . . . , a„ are arbitrary, we see that (j2ip is equivalent to (jl7p and thus to (jl6p 

We next prove the uniqueness of the solutions to the Schwinger-Dyson equation in a small 
parameters regime. 

Lemma 2. For any c > and K > i{c), there exists t{c, K) so that when \ti\ < t{c, K) for all 
i € {1, ■ ■ ■ ,k}, there exists a unique solution tr to (jl6p which satisfies 
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tr(xx*) > 0, 

tr(l„) = 1, where !„ is the trivial path consisting of no edges starting and ending at v ^ F , 
t?((eOe)P) < KP for all e £ E+ and p € N. 



Proof. The proof follows the arguments of jGMS06[ Section 2.4] that we however repeat for 
further uses. Existence is already guaranteed. To prove uniqueness we assume we have two 
solutions Ti and T2 and let 



A{n) = sup \ti{w) - T2[w)\ 

w^L ,\w\'Cn 



with \'w\ the number of letters in the word w. Then equation (1161) implies that, if we let m 
min(Ai(s(e)Ht(e)))3 > 0, 



n-l 



A{n + 1) < 2m-i ^ A{p)K'^~^-P + A{t)A{n + D - 1) 



(22) 



with D the degree of W and if Df.oW = J2i X]j=i itj with some monomials gfj with degree 
at most — 1 , then 

A{t) = max(Ai(s(e))M(t(e)))-5 ^ \t^\kl 

i 

For 7 < 1/i^, the sum /\(7) = X!p>i 7^^(p) finite and satisfies 



^(7)< 



7 



We then choose t small enough so that 

7' 



Z\(7)+A(i)7-^+M(7) 



A(i)7 



< 1 



1 - 7L 

for some 7 g]0, which guarantees that Z\(7) = and therefore 5{n) = for all n > 0. □ 
In the next section we study Schwinger-Dyson equations for laws on V'" and P. 

2.5.5. Proof of Theorem\^ in the finite graph case. In the case i = 0, we deduce Theorem [31 
First, note that the previous lemma already entails the convergence of / M~^TT:{A^)dii^ {A), 
w G L, to the unique solution of the Schwinger-Dyson equation r. Hence, we only need to 
identify this limit. It is enough to show that if we define (Tro)ij as but with Ttq instead of 
T, it will satisfy the same Schwinger-Dyson equation. Equation (jl9[) implies, for t = 











T 






T 











p 

T 



i odd 



© 
p 
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This equation clearly has a unique solution, since it defines the maps Ti,j recursively in terms 
of Tiji with J < j'. We claim that in fact r^j = (Tro)i,j where this notation means that we 
replaced tr by Trg in the definition of Tr^ ^ . Thus (Tro)ij- is given by 




T 

where ^ TL stands for the sum of all TL diagrams. Here the i + 1-st string of P is the first 
one of the j strings connected to ^ TL. If we follow the rightmost top string of P, it will be 
connected to one of the other vertical strings of P (and, for parity reasons, it will be an odd 
string). From this we see that Tij satisfy the same recursive relation as (Tro)ij. □ 



2.6. Proof of Proposition\^ Let us now consider the case t ^ Q. We have already shown part 
(b). We have also shown that any limit point of \ /ij^'^(tr(-)) > satisfies the Schwinger- 

L ^ J M->-oo 

Dyson equation; it follows from uniqueness of the solution that all such limit points are equal. 
Thus the limit 

'tTt{X^){v) = lim ^.^'^'^ {iv{A„)) 

exists and we obtain part (a) of the Proposition. 

To verify part (c), it is sufficient to show that satisfies the Schwinger-Dyson equation: 




(23) 

To see that this equation holds, we note that it simply expresses the fact that the rightmost 
top string of P must either come back to P, or be connected to a vertex coming from a copy 
of the potential W. This shows part (c). 

Finally, for part (d), let P C T'^ be a subfactor planar algebra, and assume that W ^ V (as 
an example, one could consider TL C via the embedding B i-^ wb given by equation (jlOp ). 

We note that HQ = J^w '^wW G V and w G T, then Tt{Q){v) = trtij^ui awAn]){v) = 3t{Q){v) 
by linearity. 

The infinite power series defining ^t{Q) converges, and its finite partial sums involve values 
diagrams connecting Q to several copies of W and having no unconnected strings. Since W d P 
and P is a planar subalgebra of T'^, the values of these diagrams are the same whether they are 
evaluated in P (yielding an element of (7^)0, since no unconnected strings are present) or P^ 
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(yielding an element of {'P^)q). Since V is a, subfactor planar algebra, C = {'P)o is the subset 
of ('P^)o consisting of elements which are constant on all vertices. It follows that the value of 
^t{Q){v) does not depend on v. □ 

2.6.1. Free energy. We recap here how to get from the convergence of the tracial state that of 
the free energy 

To that end we observe that, by differentiating with respect to a, 



1=1 vev u.eLBjti) ° 



Note that Fq ^'^ goes to zero as K has been chosen large enough. Thus, by Proposition [2j for 
I til smaller than t{c, L), the last integral converges uniformly for all a G [0, 1] and therefore by 
bounded convergence theorem, we get 

1=1 II ni>Q j J' "'O 



2.6.2. The case o/ Aqo and TL for 6 > 2. The previous construction only allows a countable 
set of values of ^'s (which however contains all the possible S < 2 and accumulates at 2). To get 
all possible values of 6, we need to consider infinite graphs. However, our construction below 
relies heavily on the fact that the entries of the eigenvector n go to infinity exponentially fast 
with their distance to a distinguished vertex of the graph. We will therefore restrict ourselves to 
the graph Aqo for which we shall prove that this property holds. Since not all subfactor planar 
algebras can be embedded into the graph planar algebra of Aqo, we shall restrict ourselves to 
Temperley-Lieb algebra TL in this section. In this case, it is enough to consider the infinite 
graph Aoo since this graph possesses an eigenvalue S for any possible 5 >2 with an eigenvector 
with positive entries. We let fi be the corresponding eigenvector and denote by n e N the 
vertices of Aoo- Noting that by definition for n > 3 

/i(rt — 1) + f^{n + 1) = Sfi{n) 

we see that ii{n) « 5" as n goes to infinity. We let, as in the finite graph case, jj.^ be the law 
of the independent Ms(g) x M((e) matrices with covariance l/A/-\//i(s(e))/i(t(e)) for e G 
(which is now infinite). Edges far from the origin will have variance decreasing exponentially 
fast with the distance to the origin (recall that 6 > 1). To construct the law on the infinite 
graph, we let Sn be the sigma-algebra generated by the Ae for e € En, the set of edges with 
distance less than n from the origin. We let Vn be the vertices connected to the origin by a 
path in £"„, and L{n) (resp. Lsin)) the set of loops (resp. of Lb) with all edges in En- The 
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idea is to consider the Gibbs measure indexed by infinite graphs as the hmit of the conditional 
expectation with respect to Z'„. More precisely, we define 

d/zf ^'"(dA) = Ikl^^l^d^exp J Af^t, ^(«) E ^B,iw)TTiA^)\d^^"'idA). 
We shall prove 

Proposition 3. Let K > 2. Then, there exists t{K) > so that whenever T[iaxi<i<k \ti\ < t{K), 
for any vertex v for any loop w € L{v) there exists a limit 

tTt{Xw){v) ^ lim lim /xf^'"'^'"(tr(yl„,)) . 

n— yoo A/— >.oo 

Moreover, trt(Pg)(ti) = -S^tiB) for any Temperley-Lieb tangle B and any vertex w G V+ (resp. 
V-) if the marked point of B is in an unshaded (resp. shaded) region . Finally, 

1™ 1™ TTT^ --^logZf'^'"= V ]T2^t"^^tt'-f^ 

Proof. In fact, due again to strict convexity, the negative of the Hessian of the log-density 
of /if^'^'" is bounded below independently of n on the set ||^||oo < K. Moreover, because 
(^(s(e))^(i(e)) 2 is at least of order (5" if s(e) € Ef^, we can choose t small enough so that this 
Hessian is bounded below by the diagonal matrix which takes the value C(5" for v G En+i\En 
with some positive constant C. This is enough to guarantee Brascamp-Lieb inequalities and 
concentration of measure which do not depend on the dimension and in particular on n. In 
particular 

^f'^'"( max \\A,\\^>ie^eo)S-''^')<e-^'" 

eG-Bfc + A-Bfc 

for some £o < oo, £ > Iq and rj > 0. Moreover, we can apply concentration inequalities and 
use the fact that is a Lipschitz function of the entries with Lipschitz norm of order M~ 2 
with overwhelming probability because of the above control (since <5 > 1 we obtain absolutely 
converging sequences) to see that ([T5|) still holds under /i^ ' with a constant C{w,K) inde- 
pendent of n. Therefore, by the same arguments, we obtain the convergence of /if^'^'"(tr(A„)), 
for any loop w £ L{n) as M go to infinity. The rescaled (by fi) limit r" satisfies a Schwinger- 
Dyson equation which has a unique solution, exactly as in Lemma [5J again because the controls 
are uniform in n, the range of t for which we have uniqueness does not depend on n. Note that 
this equation is the same for r" and r™ provided we consider only loops in L[n/\m — D) where 
D is the maximal number of boundary points in the potential. Define 

A{k,p):= max \t'\w) - t"\w)\ 

We see that A{k,p) satisfies ([^ for k + p < n Am. For larger k we just bound the additional 
term uniformly by An/.,„i{t)K^^'^^''^ with 



An{t) = max ^y fi{s{e))fi{t{e)) maxy^ \ti\ki 

eg E" I ' 
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We then find that if Z\(7,p) — J2k>o j'''^{^tP)i ^'"^Y P = n A ni/2, 

^' k>m/\n/2 

This yields, if (C")~^ = 1 — mii-K-y) ~ ^(^7^^^ is positive (which is always possible if t is 
small enough for some 7 e (0,-?^^^)), 

p>mAn 

which goes to zero as n A m goes to infinity for K^5^^ < 1. This shows that tr" converges as n 
goes to infinity. The limit satisfies Schwinger-Dyson equation with appropriate S and therefore 
corresponds, when restricted to P, with The statement on the free energy is derived as in 
the previous subsection. 



3. Combinatorics of a few models 

In this section, we actually compute the generating functions of a few loop models we have just 
constructed. 



3.1. The cup matrix model. We consider the case of the matrix model with tangles i?„ = cup" 
obtained by drawing n non nested strings in a white tangle with black inside. With k finite and 
fixed, we thus consider the law corresponding to the tangles {Bi, . . . ,Bk) and the coefficients 
{ti, . . . ,tk). By Propositions [5] and [21 Bn is represented in by {J2eeE+ '^{^)^eX*)'^ , e.g. 

X,,„.„, X,,,.,,,, -^C-,") ^(».-) 



B, 



B3 = ^(,.,„) ^^^^^(^,™) , etc 



and the associated Gibbs measure is, by Example [21 

Note that the families (Ae,s(e) = v), w e V+, are independent. Since we evaluate the above 
expression at words in Z„ = (X]e s(e)=t) '''(^)^e^e) for a given w we may as well consider 

the law of Z„ for a fixed v. Recalling that has variance ^/M^(^)M^e) we find that 
Zv = J2e:s{e)=v <^{e)XeX* = Y^Y* with a M„ X (X]e:s(e)=t; ^t(e)) matrix Yy with i.i.d. centered 
Gaussian entries with covariance M~^. The law of the eigenvalues (Ai, • • • , A^(^)m) of such a 
matrix is asymptotically equivalent (since we can again by Brascamp-Lieb inequality remove 
the cutoff on A and transform it as a cutoff on Zy for some K' large enough) to 

'^/^f W-^U |A. - A, I n A?"'^'^" """e-^ ^f^<^'^^ *"^-^-) n lA.e[o,K']rfA.. 
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Since iJ2e-s{e)=v Mt{e))/My converges as M goes to infinity to S, it is classical |AGZ10[ Theorem 
2.6.1] to prove a large deviation principle for the law of the empirical measure of the Xi under 
lj^M,K jj-Qj^ which one easily derives the convergence of the empirical measure. Therefore we 
deduce that 

Lemma 3. For all K > 2, and t = {ti, . . . ,tk) small enough, all n > 0, all v G V+, 

^t{B^) = lim Aif •^(tr(P;^J) = (24) 

with Vt the only probability measure on [~K,K] which maximizes, with P{x) = X]n=i^"'^"' 

It{v) := S{v)+5 / \og\x\dv{x) + / P{x)dv{x) - / xdv{x) 



where S is the free entropy 

E{n) = / / log |x - y\dv{x)dv{y). 



We can obtain an equation for vt by writing that It{i^t) > Iti^t) with the law of a; + eh{x) 
under ft and h a smooth real- valued function. Expressing that the linear term in e must vanish 
and ultimately taking h{x) — {z — x)~^ we deduce that vt is solution of the Schwinger-Dyson 
type equation 

1 \^ f I f P'ix) - 1 
dvt{x) ] +S — rdiytix) + / di^tix) = 

— X J J x[z — X) J Z — X 

which yields with Gt{z) — J j^dvt{x),Q(z) ~ J ^ ^^IZ^ dvt{x), and c — 5 J x^^dvtix), 
Gt{zf +(^- + P'{z) - Gt{z) + Q{z) + 1 = 

so that 

Gt{z) = ^ (^(1 - P'{z)) - S - y/iziP'{z)-l) + 6)^-AzizQiz)-c)) . 

It can be shown that for t small enough, vt bas connected support and deduce a formula for Gt , 
see e.g. |BIPZ78j . Indeed, according to [AGZlOi Lemma 2.6.2], for small i^'s, I't is characterized 
by the fact that a certain strictly concave function is smaller than some constant outside of its 
support and equal to it at the boundary of the support; this is only possible if the support is 
connected. Since Gt must be analytic outside the support [a, b] of vt, the formula for Gt can be 
uniquely determined by the fact that there exists a polynomial R with degree fc — 1, and a < b 
so that 

{z{P'{z) - 1) + (5)2 - Az{zQ{z) - c) = (z - a){b - z)R{zf 
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3.2. The 0{n,m) model. Consider the potential 




considered in the planar algebra of two stitched Temperley-Lieb algebras (see ^2.1.2p and recall 
the notation of ^2.2.21 

In the finite-depth case, the random matrix model is given by considering block matrices 
indexed by pairs (e, v) where e is an edge in Fr and w is a vertex in or e is an edge in Fh 
and u is a vertex in Fr with the potential of the form 

V black vertex v red vertex 

XI ^(^)'^(/)^(e.^(/))^(/.t(e))^("e,t(/))^(/,s(e))- 

eer, /en 

li Sr = n and Si, = to, we can choose Fr (resp., /],) to be the graph with one odd vertex and 
n (resp., m) even vertices, with exactly one edge between the odd vertex and every even vertex, 
the resulting combinatorics is exactly the same as of the 0{n,m) model with the potential 

n m 

i=i j=i 

which was studied for instance in |DFGGdO| . in connection with the problem of enumerating 
meanders. 



3.3. The double cup matrix model. The potential is a linear combination of two tangles; the 
tangle with two cups with black inside (with coefficient A) and the tangle with two cups with 
white inside (with coefficient B respectively). We denote by V the element of the planar algebra 
associated to it, namely (see Example [5]): 



Diagrammatically, this corresponds to tangles of the form 



where the only difference between the two pictures is that on the left, v Cz Vj,- and u,w e VL, 
and vice versa on the right. 

The associated Gibbs measure fi^^'^ (dA) is given by 
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To analyze the asymptotics of this measure, we shaU first perform a Hubbard-Stratonovitch 
transformation and then study the resuhing Gibbs measure. We first relate this measure with 
our problem. For the sake of simplicity, we restrict ourselves to (S's corresponding to finite 
graphs, and even further to the graphs A„ with n vertices and n — 1 edges. This is enough to 
characterize the generating functions by analyticity (since the index 5 corresponding to these 
graphs has an accumulation point at 2). In fact the restriction to finite graphs allows us to avoid 
dealing with a Gibbs measure on infinitely many matrices whereas the restriction to A„ ensures 
the uniqueness of the minimizers of the entropy described in Proposition [SJ and therefore their 
characterization. 

3.3.1. The partition Junction and an auxiliary matrix model. Let us first consider the partition 
function 

and introduce independent x matrices from the GUE with covariance l/M^. Then, 
assuming tA,tB positive and putting a{v) = ^/8t2 if u G V-, a{v) = \/Sti if w € V+, — 
n{v)y^ Mv{^{v)M)^^ (which approximately equals /i(i')), we get 

^M,i^ ^ J l||^||^^.e-fi:„ev"(-)MW'Tr(G. -(e)A.A:)^M(^^^^^)_ 

Note that again for a{v) small enough, the integral has a strictly log-concave density and 
therefore Brascamp-Lieb inequalities show that the matrices G„ are also bounded with large 
probability and so there exists K' large enough so that, we have 

where we used the standard notation Am ~ Bm, for two sequences Am,Bm, for a shorthand 
for AmB^ converges to one as M goes to infinity. Diagrammatically this corresponds to the 
following "breaking" of the tangles: 

We next integrate over the matrices X^^ , recalling that the entries of Ae have covariance 
(Ms(g-)Mt(g'))~^/^. Up to a constant, this provides the term 

J~J g-Tr(8iTr(log(/+a(s(e))'/(8iG,(e)+a(t(e))'Gt(e)«i/)) 

where we noticed that / e~'^^ dx = \/27r7^2 and the matrices have complex entries (so that 
each term appears twice). Here 




a(.(e))'a(s(e))-i =a(t(e))'a(t(e))-i = {pit{e)y / ^^{t{e))y^' 
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is approximately equal to one. We can finally diagonalize the matrices Gy to get 

Ml, 



7M,K 



n l|A»|<K'Z\(A'^)dA"exp 



. v^V i=l 



- log(l + a(s(e))Af'+a(t(e))Af)) 



l<«<A-fs(e) 
l<i<J\ft(,) 



with A" = (A^, . . . , XIjJ the eigenvalues of Gy and A{X') = l\i<i^j<M„ \K - I- The asymp- 

totics of log Z( ' can be obtained via the global asymptotics of the eigenvalues (At,, v G V), 
that is the convergence 



v+ 



under the associated Gibbs measure P*^'^(c?A) given by 
TT Z\(A")dA" • TT exp - 



^^n^(^")^^"- n -EE ^og{l + a{sie))xf^ 



■ exp 



A/„ 



In the sequel, we shall prove not only this convergence but the existence of two probability 
measures and vj^ so that 

xPdi^y{x) ^ J xPdiy±{x) peN, veV±. (25) 

Before attacking this question, let us summarize what information the auxiliary probability 
measure P^^'^ and (P5|) tells us about our original question. Recall that ^ is the tracial states 
constructed with the two tangles with two strings and opposite shading. We claim that i'± gives 
the law of a cup under in the following sense. 

Proposition 4. Assume (|25p and recall that a — \/8t2. Let Bn be the tangle with n non nested 
strings and black shading inside and put for small 7 



C(7,AB)-E^"-^(^") 



ri>0 



and M{z) = / X]n>o •^"•^"'^^+(^)- Then, 7(2:) = j^J^a j^^^-, is invertible from a neighborhood of 
the origin into a neighborhood of the origin, with inverse 2(7) and 



C(7,AS)-^Af(z(7)) " 



72(7) 
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Proof. Observe first that P^ ' is tlie law of the eigenvalues of {G^ ,v G V) under 

z.f'^(dAdG) = ta^exp| -^^a(t;)tr(G, ^ a{e)A,,AlM {dA,dG). 

I v(^V s(e)=v J 

Adapting the previous considerations we see that G„ J2e:s(e)=v '^{^)^eX* correspond to an 
element of a planar algebra with one string and one strip, the strip being in the white shading 
iff V G both being independent in the sense that they can not be ghicd together, and the 
strips requiring to be coupled with another strip corresponding to the same vertex, the weight 
being one. We can also show the convergence (in a small parameters regime) of the law fj^'^ 
restricted to the planar algebra generated by elements with non-crossing strings and strips. 
We denote the limit by trj. ^ corresponds to the case where we restrict ourselves to elements 
with only strings (since then expectation over the strips, that is the Gaussian variables can 
be taken) whereas i'± corresponds to restricting ourselves to element with strips only (inside a 
white or a black shading). To relate both states let us consider the expectation of an element 
Bn,p with n non nested cups with black shading inside, followed by p strips in the white region. 
Let C{p, n, £, k) be the number of possible configurations build above this tangle with I (resp. 
k) tangles with one string and one strip in the white (resp. black) shaded region. We get an 
induction relation by gluing the first strip which yields for p > 1, 

P-2 

C(p,n,£,fc)= J2 J2^k"Ce'C{pi,0,h,ki)Cip-pi-2,n,£-h,k-ki) 
pi=oei<e ki<k 

+eC{p-l,n + l,e-l,k). 

























■ 


n 


V 



The first term appears when the strip is glued with another strip of the tangle whereas the 
second one shows up when it is glued with a strip of an element with a strip and a string with 
the right shading. 
We let 

C{z, 7, a,li) = Y, C{p, n, £, k) 



l\k\ 



and conclude from the induction relation that 

~2 



az , 



C{z, 7, a, p) - G(0, 7, a, p) = z'C{z, 0, a, p)C{z, 7, a, /3) + —{C{z, 7, a, /3) - C{z, 0, a, p)) 

7 



which gives: 



C{z,j,a,p) 



7(7(0, 7, a, /3) — azC(z, 0, a, /3) 
-y — az — 'yz'^C{z, 0, a, /3) 
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Since C{z, 7, a, j3) is analytic in z, 7 small enough, we deduce that if the denominator vanishes 
so that 

7 = (1 - z^C{z, 0, a, (3))-^az =: 7(2) 

then the numerator must vanish too. Therefore, with 2(7) the inverse of 7(z) (which exists by 
the implicit function theorem in a neighborhood of the origin) we deduce 

C(0,7,a,/3) = ^C(z(7),0,a,/3). 

7 

Since if we choose a = y/8t2, (3 — y/8ti, we have C(0, 7, a, f3) = C{j, A, B) and C(z, 0, a, (i) = 
M(z), we have proved the proposition. 



3-4-. Solving the auxiliary matrix model. In this section we study the law p*^'^ and prove (j25l) 



3.4.1. Large deviation estimates and limit points. Using standard large deviation theory (see 
[AGZlOi Section 2.6]), and putting 

v):^ II log \x - y\dfi{x)dv{y), S{fi) := S{fi, fi) , 



we deduce that 

Proposition 5. Set a{v) = \/8t2 (resp. — y/8ti) if v E V- (resp. v eV+). 
and j3 = ^8ti. The law of the spectral measures 



of the matrices G„,w G V under P^ ^'^ satisfies a large deviation principle in the scale with 
good rate function 

Itivv.v G ^) ^ X! I x^dv^{x) ~ ^ ^Ji{vfS{v^) 

vev •' vev 

+ X] /^(^(^))A'('5(e)) / + a{s{e))x + a{t{e))y)dvs{e){x)dvt[e){y)- 



eeE^ 



Take F = A„,ri > 1. Then, It achieves its minimal value at a unique point so that Vy — v± if 
V E V± with {v^, v^) the unique minimizer of 



^^(^ j x^dv^{x) - S{v^ 
In particular, for any p>Q, we have 

M, 



)]+d // log\l + ax + I3y\diy+{x)diy^{y). 



lim E 

M->-oo 



^E(Arr = / x^di^^ix) tfvev± 
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Remark 1. By large deviation teelmiques as in jBAG97) . it is easy to check that {vj^.v^) is also 
the limit of the spectral measures of the two Hermitian M x M random matrices G+ , G_ with 
joint law given by 



exp [STr (g) Tr log(/ + aI(g>G++ PG^ ® I)] exp 



_^Tr(GUG^-) 



dG+dG- 



when K' is large enough. This last formula is easily obtained for integer values of 5 starting 
from the graphical rules of the double cup matrix model, and usually in the physics literature 
such expressions are then analytically continued to 5 non integer. 

Proof of Proposition [5j The first point is to show the uniqueness of the minimizers of /(. 
We let Vy be the law of —ax (resp. 1 + (5y) under for u G (resp. v e V-). Then, we have 
to minimize 

It{i^v,ve V) := H{Dy,ve V)+L{Dy,ve V) 
with, for probability measures py,v V on the real line, 

and 

L{p.,ve V):=\Y. {^f [ x'dp^ix) + Y: (^)' / (1 - xfdp^ix). 

Since /( is a good rate function, it has compact level sets (see the case V = {0} in [BAG97J) 
and therefore It achieves its maximal value. We next prove that its maximizer is unique. Note 
that L is linear in the measures. We shall prove that H is strictly convex. Indeed, put 

d{v) = tt{e e E+ : s{e) = v} + tJ{e G E+ : t{e) = v} 

and observe that when F C Aoo the degree d{v) of each vertex is bounded by one (for the 
boundary points) or by two. Therefore, the quadratic form 

veV eG-E+ e£E+ v£V 

is positive definite. We let {"jijVi) be the eigenvalues and eigenvectors of the corresponding 
matrix in R+* x RI^I and write 

Vi{u)^{u)pu)- 

Finally S is strictly log-concave on the set of real valued measures with fixed mass, as was 
proved in |BAG97| . Applying it to the measure pi — '^^^^y ^J'(u)vi{u)pu with given mass 
^i^gy Wi(u)/x(u) and using that the 7i's are positive, we deduce that H is strictly convex. 
Therefore, It (and by the same argument S) achieves its minimal value at a unique point 
[vy^v e V) (resp. {vj^,v-)). We next show that it has to be Vy = v± \i v ^ V±. In fact, the 
infimum of It is characterized by the fact that for all v g , the function /„ (x) given by 



2^(w) j \og\x ~ y\dvy{y) - ^ A*(^(e)) j \og{l + ax + I3y)dvt{e){y) 



x^ 



e:s{e)—v 
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is constant on the support of i^y and non positive outside. We have the same equation for 
V G E- with a and /3 exchanged. The same characterization holds for which is such that 
the function /_|_ given by 

-S J log{l + ax + I3y)diy-{y) + 2 j \og\x - y\dv+{y) - ^ 

is constant on the support of and non positive outside (and again a similar equation for 
j/_ with a and /3 exchanged). Putting v.^ = h'± for v (z V±, we find that fy = tJ-{v)f± as 
J2e-s(e)=v A'(^) ~ '^/^(^)j 9,nd therefore is indeed constant on the support of and non positive 
outside. Thus Vy = v± for v & V± \s & solution, and by the previous argument the unique 
solution. □ 



3.4-2. Properties of the minimizers of the rate function. We can give the following characteri- 
zation of the minimizer (i/-|_, j/_) of S. 

Lemma 4. Let a — \/8tA and (3 — ^/StB. There exists to > so that for \t\ < to, 

— v± has a connected support S± included in [—2 — A{t), 2 + A{t)] with A(t) going to zero as 
t goes to zero. 

— There exist functions {f±,g±) which are analytic in a neighborhood of S± and so that for 

G±{z) = J j^diy±{x) = /±(z) - v/^^ (26) 

with the branch of the square root chosen on . Moreover g± is real on the real line and 
S± = {x : g±(x) < 0}. We define G±{z + iO) and G±{z — iO) as the limit of G± when z goes 
to an element of S± from above or from below. 

— For all X (z S± we have with Q!_|_ = a, a_ = /3, 

/ l-a±a^ \ ^ ^ ^ Q^^^ ^ -Q^ + G±{x- lO). (27) 

— There exists at most one solution to (j27p so that G± are the Cauchy-Stieltjes transform of 
probability measures v± supported by [—3,3]. 

Proof. The fact that the support of and is connected is a consequence of Remark [T] 
and |GS091 Theorem 4.4 and Theorem 4.2] which asserts that the limiting measures v± have 
connected supports when the potential is strictly locally convex, which is the case when t is 
small enough (note here that the potential is not a polynomial, however it expands in absolutely 
converging power series when t is small enough and x,y are bounded by K' so that we can 
apply the techniques from |GS09| to represent the measures as the law of an element of the G* 
algebra generated by the free Brownian motion). Another way to see this is to notice that the 
function /+ is strictly concave outside of the support and continuous on the support where it 
takes only one constant value; hence the support can not be disconnected. The fact that the 
support is bounded is a direct consequence of Brascamp-Lieb inequalities. 
To deduce the second point, we obtain an equation on {v+, v-) by writing 
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with i/£ the law of a; + CfT-±{^) for bounded continuous functions h± on M. Writing that the 
hnear term in ^ must vanish results with the equation 

(28) 



x-y 
1 + ax + py 

Taking h+{x) = ~^{z + and h_{x) = (z - x^^ we get that, with m{z) = G+(-(l + 

Pz)/a), 



^m{zf + G-{zY + 5-m{z)G-{z) - zG-{z) + 1 + ^(1 - (1 + — = 



a a 



G^{z) = \{b{z)-^b{zY-Aa{z)) 



which gives 

with the cut of the square root on Mr and 

h{z) = z - 5^m{z), a{z) = ^m{zf + 1 + ^(1 - (1 + —)m{z)). 
a a a a 

But, for small t, when z is in the neighborhood of S*-, — (1 + I3z)/a is in the neighborhood of 
— 1/a which is far from the support 5*+. Hence, a{z) and b{z) are analytic in the neighborhood 
of S- and also take real values. This completes the proof of the second point. 

For the third point, it is enough to take /i+ = or /i_ = in ([25)) with the remark that 

^^^^ ' ^^^^ dv±[x)dv±[y) = / h[x)[G±[x + iO) + G±{x-iQ)]dv±{x) 
x-y J 

and use the continuity of G± above and below the cut due to the previous point to obtain the 
desired equations almost surely and then everywhere. 

For the last point, since (P7| is equivalent to (pS)). we show the uniqueness of the solution 
by taking h+{x) — {z — x)^^ to deduce that 

-l + zG±{z)-G±{zf ^5a± f -^a^G^{~^-^^^^)diy±ix) =: e±{z) 

J Z — X Qip 

SO that for z sufficiently large 

1 



G±{z)^-{z-^z^~4e±{z)) 

where we have taken the usual determination of the square root because e±(z) is small, since 
aipGip( — ^^^^) is uniformly close to one for x G S± and t small. Therefore, if we have two 
solutions G and G, we find that there exists M{t) finite so that for t small, 

sup \G±{z) ~ G±{z)\ < I sup \G^{z) ~ G^{z)\ 

|z|>M(t) NI>JVf(t) 

which results with G^(z) = Gip(z) for z large enough and real, and then for all z in the 
complement of S± by analyticity. 
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3.4-3. Characterization of the minimizers of the rate function. In this section we completely 
characterize the measures {v+,v-) by their Cauchy-Stieltjes transform. To simplify the no- 
tations, we let i>-4_ and v- be the law of —ax and 1 + jSx under and respectively. By 
Lemma m for t small enough, v+ and j>_ have disjoint connected supports [01,02] and [61, ^2] 
around the origin and the unity respectively. Our study will be restricted to this case, which 
therefore include small t's but eventually a larger class of parameters. We will first proceed by a 
reparametrization of the Cauchy-Stieltjes transforms of Z/+ and , which will allow to obtain 
very simple equations, and then solve these equations. 

• Parameterization. Consider the standard parameterization of the elliptic curve = 
(z - ai){z - a2){z - bi)(z - 62): 

uiz)^-^ih-a,){b,-a,) / , .. , ,u, (^9) 

^ Jb2 y/(z' - ai){z' - a2)[z' - bi)[z' - b2) 

where the path of integration avoids the segment [01,62] (for z e [02,61], we choose to come 
from the upper half plane, though this choice is irrelevant, see the comment below on 2K- 
periodicity), and the square root in the denominator is defined as having cuts [01,02] and 
[61, 62] and such that at infinity it behaves like z^. 

The image of z G C U {00} — ([01,02] U [61,62]) u{z) = {Re u{z),lmu{z)) is a rectangle 
of the form [—K, K[x]0, iK'[, where K and K' are usually called quarter-periods (though they 
will be half-periods in what follows). Now, note that crossing the line [02,61] corresponds to 
u — > U + 2K. Since all the functions of z we shall consider are smooth when one crosses this line, 
they can be made into periodic functions on the strip < Imu < K' . The map z 1— >■ u{z) is then 
an analytic isomorphism from the Riemann sphere C U {00} minus the two segments [01,02] 
and [61,62] to ]K/2KZx]0,iK'[. If one extends the map to the cuts [01,02] and [61,62], then 
the result depends on whether one approaches them from the upper or lower half-plane, and 
they get sent to Imu = K' et Imu = as described on the figure. More precisely, if o €]oi, 02[, 
then u(o + iO) + u{a — iO) — 2iK' (the two images of o are symmetric w.r.t. iK' on the line 
Imu = K'), and similarly if 6 G]6i, 62[, then u(6 + iO) + u(6 — iO) = 0. 

iK' 




The inverse map z(u) can be expressed in terms of Jacobi's elliptic function sn, and can be 
deduced from the following identity: 



2/ i2\ oi - 61 62 - z 



sn^(u,fc'') = 

62 — 61 Ol — z 

with e - 



)l { ■ Note finally that 

(62 — a2)(oi— ai) -J 



y . . * f^°° V(6l-Oi)(62-02)rfx 

lim u[z) = - — —: Uoo G 

2 Jb^ V(a; - oi)(a; - 02)(a; - 6i)(x - 62) 
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is a pole of z{u). It is a simple pole as one easily sees that for z large u{z) ^ Uoo + + o{l/ z) 
with a non vanishing constant z_i. Moreover, z{u) is analytic everywhere else in [—K\K] x 
[0, iK'] by the implicit function theorem, and, once analytically continued to the whole complex 
plane, is even and elliptic (doubly periodic), with poles at ±Uoo (mod 2K,2iK'). 

• Resolvents. We let uj± be the reparametrization of the Cauchy-Stieltjes transform of 
and V- respectively; 

^±(w) = / -TT diy±(x) 

J z[u)-x 

The functions uj±{u) are analytic in the strip < Imu < iK', and according to the second 
point of Lemma m they can in fact be analytically continued to some neighborhood of the 
closed strip < Imu < iK' . Indeed, if < | Rezi(z)| < K the mapping z H' u{z) is invertible 
around z and the extension in the z variable translates directly into the u variable. If z = oi, 02 
(resp. 61,62 for G_), one has according to (with a more detailed study of g± which shows 
that for t small enough g'^ does not vanish in a neighborhood of S±), G±{z) oc {z — z')^/^ as 
z' — >■ z; but this matches the behavior of u{z) in ([2^ . so once again Ld±(u) is a well-defined 
analytic function in the neighborhood of 0, iK', ±K, iK' ± K. 
By (P7)) . we have 

u}+{u)+uj+{2iK' -u) - 5u}^{u) = P+{u) lmu = K' (30) 

u)^(u) + uj^{—u) — 5uj-f-(u) = P^(u) ImM = (31) 

where P+{u) = z(u)/a and P-(u) — {z{u) — 1)//?. P± are even elliptic functions (with periods 
2K, 2iK'). 

We also have the following additional conditions: 

W-i-(u) — a;+(— m) lmu = (32) 

a;_ {u) = w_ {2iK' - u) Im u ^ K' (33) 

expressing the fact that the Cauchy-Stieltjes transform of v+ is analytic in a neighborhood of 
[61, 62], so its values at z ± iO, z e [61, 62] should be equal; and similarly for v^. 

Now these equations can be repeatedly used to extend u)± to the whole complex plane: for 
example u — ^ 2iK' — u maps the strip < Imu < K' to the strip K' < Imu < 2K' , so we 
can use Eq. (j30[) as a definition of ujj^ in this new strip and equation (j30l) precisely ensure that 
the two definitions coincide at their common boundary Imu = K'; and so on. This way we 
obtain meromorphic functions w±(u) defined on the whole complex plane, and by uniqueness 
of analytic functions that coincide on a set with accumulation points, we deduce that the above 
equations are true for all u: 



uj+{u + 2K) ^ uj+{u) (34) 

Lo^{u + 2K) ^ uj^{u) (35) 

a;+(u) = a;+(— u) (36) 

uj^{u) ^ uj^{2iK' -u) (37) 

UJ+ (u) + UJ+ {2iK' - u) - Suj- (u) = P+ (u) (38) 

w_(u) + u;_(-u) - 5a;+(u) = P_(u) (39) 



Furthermore, these functions must possess a zero at Uqo '■— u{z = 00) and a prescribed derivative 
at Uoo- 
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• General solution of the saddle point equations. Now we have a well-posed analytic 
problem, which can be solved explicitly. Set S = q + q'^ . li \6\ <2,qis not real and has modulus 
one, e.g. if J = 2cos(7r/n), q = e*'^/". If 5 > 2, 5 is real and can be chosen in ]0, 1[. Set 

(p±{u) — q^^uj+{u) — — R±{u) 

where R±{u) — {q^^ P-i-{u) + q^'^P^{u)). Then the equations can be recombined into: 

ip±{u + 2K) ^ip±{u) (40) 
(p±{u + 2tK'):^q'^'^ip±{u) (41) 

where the first point is a direct consequence of (IMl) and whereas the second is obtained by 
multiplying psp by q^ and (j39p by q^^ and adding the two corresponding equations. Moreover, 
((55)) and ^ imphes that 

(p±{~u) ^ -^^{u). (42) 

Thus we may consider (p+ only. Furthermore we know that the only poles of ip± in the fun- 
damental domain [—K; K] x [—iK'; iK'] are at ±Moo i.e. z — >■ 00; they appear because of the 
inhomogeneous terms i?±, which have such poles. We can therefore express 1^3+ (it) in terms of 
9 functions. Define to be the rescaled 9i function, or explicitly 

0{u) = 2 f2 e-("+i/2)^-i^7i^ sin(2n + 1)^ (43) 



2K 

n=0 



which satisfies 



0{u + 2K) = -0{u) 
0{u + 2iK') = -e^(^'-™^/-^6>(w) 

and with a unique simple zero at u = (mod 2K^ 2iK'). 
Then we have 

Proposition 6. Write q — e*'^" with v real if S < 2 and purely imaginary if S > 2. Then, 

. . e{u~u^-2vK) 0{u + Uoo~2vK) 
(p+(u) ^ c+ — 7 h c_ — 7 (44) 

^ ^ O{u-Uoo) O{u + Uoo) 

where 

if z{u) ~ z^i/{u — Woo) + 0{1) as u goes to Uao- 
cannot 

Proof. This can be viewed as a consequence of the Riemann-Roch theorem, but we give here 
an elementary proof. 

If J > 2, let us first rule out the possibility that q^ = e"^'^"^ for some integer number 
n. Indeed then e~^™^/^ (p+{u) is also elliptic and therefore the sum of its residues vanishes. 
But the latter is given by e-™'^"-/-^(Q!-i + qp-^) + e™^«~/^(a-i + q^^l3^^) which cannot 
vanish. Thus, 7^ g-27mK /k ^^^^ therefore 0{2vK) 7^ (which allows in particular to define 
c± above). 
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Next, let denote the right hand side of (HH) and observe that by the properties of the 
functions 0, if>+ satisfies (|4T|) and (|40|. 

The formula for c± is obtained by requiring that (p^ have the same residues at u ~ ±Moo as 
our solution tp^. Indeed, 0{u) ^ 0'{Q)u as u goes to zero so that shows that 

~ , N 0{2vK) ^, , 

ip+(u) c+ + O 1 

0'{{))[U - Uoo) 

whereas both lo^ and a;_ go to zero and 

R+{u) - —^{a-^q + (i-^q^)z{u) 



1 - g2 z^i{l - q'^){u - Uoo)' 

The formula iy9+ — quj^ — — i?+ allows to conclude. The same reasoning works as u — — Uoo 
by using = -(p-{u). 

Let us finally show that / := tp+ — (p+ must vanish. Indeed / is holomorphic and therefore 
g := f'/f is holomorphic except where / vanishes, where it has only simple poles, with non- 
negative residues. But since / satisfies (HT|) and (HUl) . g is elliptic and therefore the sum of its 
residues vanishes. Hence, the residues of g vanish, and therefore by Liouville's Theorem, g is 
constant, resulting with f{u) = e^" for some constant 7. But then implies that 7 = inn/K 
and q^ = e^^TrriK /K ^ which we excluded earlier. 

Finally, we need to fix the parameters oi, 02, 61, &2- 

The first way is to notice that G±{z) is an analytic function in a,P,z in a neighborhood 
of the origin as, by Remark [U it is the Stieltjes function of the limiting spectral measure of a 
matrix model with strictly log-concave density which, even though not polynomial, expands as 
a power series, see [GS09| . The coefficients of these series can be computed recursively by the 
Schwinger-Dyson equation. Finally, by (pS)) . the boundary of the support [ai, are determined 
by g+{ai) = 0, i = 1, 2 which shows that there is a polynomial P so that P{ai, G+{ai), a, f3) — 0. 
The implicit function theorem then implies that is an analytic function of a,/3 for i — 1,2 
whose expansion can be deduced from the expansion of G^. The same applies for 61, 62- 

The second way to determine these boundary points uses the explicit formula in terms of 9 
functions and the reparametrization of the problem in terms of p :— exp(—TTK'/K) and of 

Note that because (ai, a2, &i, 62) expand analytically in a, (3, so do {K, K' , Uao) with 

fm>-l 

K = u{h) - ^",™/?"""' and K' = -m(ai) = J] 



^00 — 

n+m>— 1 



As a consequence, (p, k) also expand in terms of (a,/?), with k ^ yoJP and p ^ \foip when 
a, /3 are small but a/ P oi order one. Again by the implicit function theorem, we can invert this 
expansion and obtain a, /3 as a power series in (p, k), and therefore also (ai, 02, foi, 62), Z-i and 
ziu). We can then identify the expansion of {q — q^^)u)+ and + -R+ — (p-R- around Uqo 
to compute w+ recursively. In appendix A, we provide the first few orders of the power series 
expansion of some quantities as a function of p and n, and their diagrammatic meaning. 

Proposition 7. If 6 = 2cos(7r/7i), n > 3, G± satisfy an algebraic equation. 
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Proof. Observe that since q = e*'^/", by equations (j40l4ip . the functions ip± are elhptic with 
periods {2K,2niK'), and therefore so are uj±. The function z{u) is also ehiptic with these 
periods. But a fundamental theorem of elliptic functions WW96, section 20.54] states that two 
elliptic functions with the same periods are related by an algebraic equation: there exist two 
polynomials P± so that 

Composing with u{z) shows the existence of an algebraic relation. 

We can in fact determine the degree of P± by a slightly more explicit construction of these 
polynomials; we find it is at most 2n — 2 in z and 2n in uj. 

This result is to be compared with Theorem 15 of jBBMOQj . Indeed, their generating series 
M[q,v,t,w, z]x,y) is closely related to our G±{z;a, I3;5). The correspondence of parameters 
goes as follows: q = ^ v — 1 -\- 5 a / ji\ among the three parameters t, w, z, one is redundant and 
if z is set to 1 one has t = /S, but w is fixed to be w = 1/S {w and z are parameters weighing in 
our language white and black regions between tangles; note that introducing an extra parameter 
in our matrix model to let w vary is possible and would make no difference in the exact solution 
presented so far, so that Prop. [7] would still hold). x,y are "boundary" parameters similar to 
our parameter z, in the sense that they give a weight to a particular edge (or vertex, or face) 
of the planar map. However they are not exactly the same, and therefore direct identification 
of our generating series is not possible; only specific identities can be written, such as 

- / xdD+ix) ^ M{S^, 1 + Sa/P, 13, 1/S, 
a J 

where the left hand side is the l/z^ term in the z — > cx) expansion of G+(z), and the right hand 
side is the generating series of |BBM09| with certain specializations of its parameters. 

A. Analyticity of J?' as a function of fugacity 

Let Si, . . . , Sk and S be elements of the Temperley-Licb planar algebra with fugacity S and set 
for complex parameters ti, . . . ,tk, 

■^tAS)= E n^^i E <5#'°°P-"^ (45) 

ni,...,nfe=0 1=1 PeP{ni,...,nk,S) 

where we sum over all admissible planar maps built on Si, . . . , Sk, S . Then we state that if hi 
is the number of boundary points of Si, then 

Lemma 5. There exists a positive constant B so that for all ti, . . . ,tk,5 G C'''^^ so that 
maxi<i<fe |(5|"2"|t,j| < B, ,%_s{S) is a well defined absolutely converging series, and therefore 
t,S ^ ^t,s{S) is analytic on this set. 

Indeed, the number of loops is bounded by one half of the total number of boundary points, i.e., 
i(6 + ^ nibi), whereas before hi is the number of boundary points of Si and h is the number 
of boundary points of S. Therefore, the coefficients of the series are simply bounded by 

6(ni,...,nfc) := j| ■ 

i=i 
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and the sum can be enlarged to all planar maps that can be built over (resp. one) vertices 
with degree given by the number of boundary points of Si, 1 < i < k (resp. S). It is well known, 
see e.g. |GMS06. p. 255], that the number of such maps grows as ni=i ^^J^"^"' for some 

finite constant A. Hence, ,%^s{S) is an absolutely converging series on maxi<i<fc < 1, 

domain on which it is analytic. 

Note that it is expected that as one increases the ti, one should eventually reach a hyper- 
surface of singularities which signals the boundary of the analyticity region in the variables 
t — (ti,. . . ,tk)- This singularity is usually present in matrix models and is explained by the 
proliferation of planar maps: typically the number of planar maps grows exponentially with 
its number of vertices and this produces a finite radius of convergence of the corresponding 
generating series. This is certainly what happens in the cases studied in section[31 The model of 
13. H is closely related to the so-called one-matrix-model, whose possible "critical behaviors" (i.e., 
types of singularities) are well-known |DFGZJ95j . No exact solution is known for the model of 
13.21 but a conjecture on its critical exponent is proposed in [DFGGOO] . Finally, it is expected 
that the model of 13. 31 has a critical behavior of the type of the 0{n) model on the line A = B, 
and that of pure gravity on other lines of constant ratio A/B ^ would be interesting to find 
an interpretation of these critical behaviors in the present context, i.e. in terms of properties 
of ^t. 
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B. First few diagrams of the double cup matrix model 



The parameters ai, 02, 61, 62 of the auxiliary model of sect. 13.41 can, as already mentioned, 
be fixed by appropriate expansion of uj+{u) around ±1*00- In practice, the resulting equations 
reduce to first degree equations on the condition that one solves them parametrically in terms 
of the elliptic nome p := exp^-^TrK' / K) and of Uoo- All other quantities can then be obtained 
in the same parameterization. 

As a check, we shall here write the first few orders of the expansion for small a, j3. The 
correct scaling as a, /3 at fixed ratio is to keep the quantity 

fixed while sending p to zero, so that a//3 ^ and a/3 ^ p^ . We find the following expansions: 
a = K p- (k(35 + 2) + K"i(2^ + Q))p^ 

+ ((8(52 + 55 _^ 3)^2 _^ (^g^2 ^ ^ 24) + (5(5^ + 125 + IT)^"^) p^ + 0(/) 
P ^ \p - (k-1(3(5 + 2) + k{25 + &))p^ 

+ ((8(5^ + 5(5 + 3)k"2 + (8(5^ + 45(5 + 24) + (5^^ + 12(5 + l?)^^) p^ + ©(p"^) 



02 = K 



-2pl/2 + + 2((3 + 5)k-^ + (1 + (5)^2)^3/2 ^ 0(p2) 
2pl/2 + ,5p - 2((3 + (5)Ai-2 + (1 + ^)K2y/2 ^ Q(p2) 



61 = 1 — K 



62 = 1- 



2pl/2 + ,5p - 2((f + (5)^-2 + (3 + 5)^2)^3/2 + (9(p2) 
-2//2 + + 2((1 + 5)k-2 + (3 + 5)^:2)^3/2 ^ 0(^2) 

xdv+{x) = n 5p-5{n{25 + 1) + k"^(5 + 5))^^ 
+ (5 (k2(452 3j2 + 24(5 + 1 + K-2(-2^2 + 4^ _^ ^3 ^ (^(-^4) 

a;2 dD+{x) = k\p+{k{5'^ - 2(5 - 2) - 2k-1(3 + (5))^^ 

+ {k^{-A5^ + 3^2 + 35 + 3) + (-2(5-'^ - 4(5^ + 36(5 + 24) 
+ «-r2(5(52 + 12(5 + n))p^ + 0(p*) 

(the expansions of ai, a2, 61, 62 are only given up to order p3/2 because of issues of space). 
Inverting the first two expansions and inserting the result in the last two yields 



xdv+{x) = (5 + (5(1 + 5){a + fi) 
a J 

+ (5((2 + 55 + 25'^)a^ + (6 + 8(5 + 4(5^)^/3 + (2 + 5(5 + 2(52)/32) + 

Similarly, 

-^^y x2dz:^+(x) - a) =5(l + (5) + (5(2 + 55 + 2(52)a + (5(3 + 4(5 + 252)/3 + ... 
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+ a5 f^^l +135 
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xdv+{x) = 5 

a J 




2x2 



1 



K3VVQ 



2x3 
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Similarly, we have: 
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C. The stitching of two plcinar algebras. 

Let V = and Q = be two subfactor planar algebras (assumed spherical for simplicity). 
We define a new subfactor planar algebra 

V®Q = {P®Q)t 

which will not be irreducible even if both V and Q are. 

The vector spaces oiP @ Q are defined as follows. Fix the even number 2n and consider the 
possible partititions tt of {1, 2, 3, 2n} into two subsets of even sizes 2(p^) and 2{qTr), whose 
elements we call of type V and Q respectively. Then 

(P ®Q)t = ^iP® Q)t with (p © Q)t = P± ® Qt ■ 

TT 

To define the action of planar tangles on P © Q we use the notion of a string-coloring. Given 
a planar tangle T with discs Di as in section 2.1, a string-colouring cr of T is an assignment 
of P or Q to every string of T so that if one removes all the strings of either color, one gets 
two planar tangles T-p and Tq when one takes as initial segments, with shading, the intervals 
containing the initial segments of T. In particular there must be an even number of strings of 
each color incident to every disc of T, but this is not a sufficient condition for the shadings to 
be coherent. 

A string-coloring cr of T defines: 

a. a partition Wa of the boundary points for each disc (numbered from 1 to 2bj starting at the 

first boimdary point after the initial segment in clockwise ordcr)into V and Q points, and 

b. two planar tangles T-p and Tq by removing all the strings of the other color and taking as 
initial segments for discs the ones containing the initial segments of T. The shadings of T-p 
and Tq are determined by that of the initial segment of the outside boundary of T. 

By multilincarity it suffices to define the action Zt of a planar tangle T, with k internal 
discs, on a fc— tuple {xi, X2, Xk) of elements oiP ® Q where Xi = (g) with Vi G 

and WiQf^. 

Suppose such an element Xi of {P © Q)^, is assigned to each Di, then we define 

Mt (xi , X2 , . . . , Xfc ) = ^ Mt^ {v^- ,v;\...,vl')(g, Mt^ ,w^',...,wl'') 

a 

where a runs over all the string-colorings of T. 

It is clear that this action is compatible with the gluings. 

The ^-structure on P © Q is derived in the obvious way from those of V and Q. 

Notes, (i) Another way to describe the action is as follows: suppose elements Vi(^Wi € (P © Q)^. 
are assigned to the internal discs Di of T. Then the value of Mt is zero unless the colouring of 
the boundary points implied by the tt^ extends to a string colouring of T and then this value 
is the sum over all such extensions a of Mt^{vi,V2, ■■■Vk) ® Mt^{w\,W2, Wfe). If there are 
strings connecting the outside boundary of T to itself this element will lie in more than one 
direct summand of (P © Q)^- 

(ii) It is clear from (i) that the loop parameter of P © Q is the sum of the loop parameters 
of V and Q respectively. 
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(iii) Positive definiteness of the inner product is also clear from (i)-tlie various (P © Q)^ 
are orthogonal and in each one the inner product is just the tensor product inner product for 

(iv) It is the exponential generating functions for V and Q that behaves well under this 
operation. 

(v) For the inductive limit algebra structure ofV ® Q a, complete set of centrally orthogonal 
minimal projections is given by the tensor products of such sets of projections for V and Q. 
Thus the vertices of the principal graph of T' © Q is the Cartesian prodiict of the vertices of 
the principal graphs of V and Q and {p, q) is adjacent to {p' , q') iff p is adjacent to p' or q is 
adjacent to q' but not both. For instance if P has principal graph A3 and Q has principal graph 
A4, then that ofV @ Q is: 



A3 



A4 

(vi) We see that in the loop basis description off "P © Q, from each vertex of a loop one may 
choose to travel on the principal graph of V or that of Q. 

(vii) If we were dealing with unshaded planar algebras we would simply remove the restric- 
tions on the parity of the numbers of boundary points and the partitions tt. 

(viii) A TL basis diagram in V ® Q consists of a sum over all planar partitions of the 
boundary points into V points and Q points with that basis diagram regarded as a tensor 
product of its V part and its Q part. 
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